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1.  SUMMARY 


A  novel  variational  method  to  evaluate  the  stiffness  reduction  and  the  internal  stresses  in 

[0°  /90°  laminates  has  been  developed  and  has  been  applied  to  analysis  of  glass/epoxy  and 
m  d  is 

graphite/epoxy  laminates.  The  method  has  been  extended  to  evaluate  thermal  expaiy&on  coefficients 
of  cracked  laminates  and  also  internal  stresses  due  to  temperature  change.  Results  demonstrate  the 
stiffness  .Poissons  ratio  and  thermal  expansion  coefficients  changes  with  crack  density.  In  particular, 
it  has  been  shown  that  the  values  of  thermal  expansion  coefficients  of  cracked  laminates  also  depend 
on  the  signs  of  load  and  temperature  change.  Internal  stresses  obtained  convey  important  information 
about  sources  of  continued  internal  failure. 

A  comprehensive  survey  on  damage  in  fiber  composite  materials  has  been  prepared  including 
classification  of  the  various  kinds  of  damage  and  description  of  methods  of  analysis. 

The  differential  scheme  approximation  for  effective  properties  of  composite  materials  has  been 
modified  and  generalized  to  the  case  of  cracked  materials.  The  method  is  of  general  nature-Specific 
results  have  been  given  for  stiffness  reduction  due  to  cracks  for  isotropic  materials  containing  many 
randomly  oriented  elliptical  or  penny  shaped  cracks  and  for  orthotropic  (fiber  composite)  sheets 
containing  many  aligned  cracks  (along  fibers).  The  method  corrects  serious  shortcomings  of  the  usu¬ 
ally  employed  self  consistent  scheme. 
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2.  RESEARCH  OBJECTIVES 


(a)  To  develop  analysis  of  stiffness  reduction  of  laminates  due  to  intralaminar  crack  distribution 
in  more  than  one  direction. 


(b)  To  develop  analysis  of  change  of  thermal  expansion  coefficients  of  laminates  due  to  intra- 
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laminar  crack  distributions. 

(c)  To  evaluate  internal  stresses,  due  to  load  and  temperature  change,  in  laminates  with  intra¬ 
laminar  crack  distributions  and  to  utilize  these  to  assess  sources  of  progressive  internal  fai¬ 
lure. 

(d)  To  evaluate  stiffness  reduction  of  solids  due  to  distributions  of  cracks  of  various  shapes. 


(e)  To  analyze  by  variational  methods  stress  concentrations  near  edges  of  laminates  and  holes. 


I 


3.  STATUS  OF  RESEARCH 


Under  a  previous  AFOSR  contract  we  have  developed  a  novel  variational  method  which  enabled 
us  to  evaluate  stiffness  reduction  and  approximate  internal  stresses  in  a  cross-ply  laminate  with  only 
one  ply  cracked.  The  stiffness  reduction  results  were  found  to  be  in  remarkable  agreement  with 
experimental  results,  [1],  Under  the  present  contract  we  have  generalized  this  method  to  cross-ply 
laminates  in  which  all  plies  may  be  cracked  along  the  fibers  42].  We  have  called  such  laminates  - 
orthogonally  cracked.  It  should  be  noted  that  to  the  best  of  our  knowledge  the  literature  contains 
only  one  publication  on  the  subject  of  analysis  of  orthogonally  cracked  laminates  which  employs 
very  extensive  numerical  analysis  while  the  present  anlysis  is  mathematically  so  simple  that  the  res¬ 
ults  can  be  evaluated  on  a  programmable  calculator. 


We  have  obtained  expressions  for  Youngs  moduli  and  Poissons  ratios  of  orthogonally  cracked 
laminates  for  arbitrary  crack  density  of  intralaminar  cracks  as  well  as  approximate  three  dimensional 
stresses  in  the  laminate.  We  have  found  that  only  the  cracks  transverse  to  the  load  have  a  significant 
effect  on  Youngs  modulus  in  that  direction  while  Poissons  ratio  is  significantly  affected  by  the  two 


crack  families.  The  internal  stresses  reveal  that  the  stress  normal  to  laminate  midplane  (peeling 
stress)  is  a  major  source  of  continued  crack  formation  ,[2]. 

To  complete  the  elastic  properties  of  the  orthogonally  cracked  laminate  we  also  need  the  shear 
modulus.  Unfortunately  our  considerable  efforts  to  evaluate  this  modulus  have  not  been  successful 
because  of  some  peculiar  mathematical  difficulties  and  we  hope  to  resolve  this  problem  at  some 
future  time. 

Other  important  properties  are  the  thermal  expansion  coefficients  of  the  cracked  laminate.  We 
were  able  to  evaluate  these  in  novel  and  simple  fashion  [3]  on  the  basis  of  some  general  theorems  on 
thermoelasticity  of  composite  materials  and  the  internal  stresses  given  in  [2],  The  results  show  that 
the  thermal  expansion  coefficients  diminish  significantly  with  increasing  crack  density  to  attain 
asymptotic  values.  Furthermore,  it  has  been  shown  that  the  values  of  the  thermal  expansion  coeffi¬ 
cients  depend  on  the  nature  of  the  load  (tension  or  compression)  and  temperature  change  (heating  or 
cooling)  since  internal  cracks  may  all  be  open  ,  some  open  some  closed,  or  all  closed. 

We  were  also  able  to  evaluate  the  internal  thermal  stresses  in  an  orthogonally  cracked  laminate 
[6]  by  utilization  of  variational  principles  of  thermoelasticity  and  admissible  stress  fields  of  the  type 
used  in  [2].  These  stresses  indicate  which  way  internal  failure  proceeds  due  to  temperature  change. 

A  very  important  problem  which  has  received  much  attention  in  the  literature  is  stiffness 
reduction  of  a  solid  by  various  distributions  of  cracks.  For  example  :  elliptical  cracks  in  different 
orientations  or  aligned  line  cracks.  Exact  solutions  for  such  problems  are  available  only  for  dilute 
concentration  of  cracks  which  implies  that  mutual  distances  between  cracks  are  so  large  that  they  do 
not  interact,  which  is  not  a  practical  situation.  The  current  most  popular  approximate  approach  to 
this  problem  is  the  self  consistent  scheme  but  the  results  of  this  method  are  often  unreliable.  We 
have  been  able  to  devise  another  approximate  method  which  is  as  general  as  the  self  consistent 
scheme  but  does  not  exhibit  its  shortcomings.  The  essence  of  the  method  is  modification  and  gener¬ 
alization  of  the  approach  known  as  the  differential  scheme  .which  has  been  employed  to  evaluate 
effective  properties  of  composite  materials,  to  cracked  materials  [5]. 
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According  to  the  contract  statement  of  work  we  have  examined  the  capability  of  our  variational 
approach  to  obtain  stress  concentrations  at  edges  of  laminates  and  holes.  After  considerable  work  we 
have  reached  the  conclusion  that  the  method  is  too  complicated  for  this  purpose  and  that  the 
appropriate  approach  to  such  problems  is  numerical  analysis  by  finite  elements. 
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ANALYSIS  OF  STIFFNESS  REDUCTION  OF  CRACKED 
CROSS-PLY  LAMINATES 

Z.  HASHIN 

Department  of  Solid  Mechanics,  Materials  and  Structures,  Tel-Aviv  University,  Tel-Aviv,  Israel 

Abstract — Stiffness  reduction  of  cracked  [0i/90j]t  laminates  is  analyzed  by  variational  methods  on 
the  basis  of  the  principle  of  minimum  complementary  energy.  For  this  purpose  admissible  stress 
systems  are  constructed  which  satisfy  equilibrium  and  all  boundary  and  interface  conditions.  The 
optimal  stress  field  is  then  determined  by  minimization  of  complementary  energy.  The  analysis  allows 
for  crack  interaction  and  random  crack  distribution.  Results  are  given  for  Young's  modulus,  shear 
modulus  and  Poisson's  ratio.  Young's  modulus  results  are  in  excellent  agreement  with  experimental 
data  for  [0°/90j],  glass/epoxy  laminate. 


1.  INTRODUCTION 

During  static  or  cyclic  loading,  fiber  composite  laminates  in  which  the  matrix  is  polymeric  develop 
distributions  of  intralaminar  cracks  which  extend  along  the  fibers,  traverse  the  plies  through  their 
thickness  and  are  essentially  perpendicular  to  the  planes  of  the  plies.  Other  kinds  of  cracks  are 
interlaminar,  i.e.  plane  cracks  within  the  interfaces  between  plies.  Such  cracks  may  be  initiated  at  the 
intersection  of  intralaminar  cracks  as  well  as  at  laminate  edges  due  to  high  interlaminar  stresses 
which  develop  at  such  locations.  Such  cracks  are  the  major  damage  which  occurs  in  a  laminate  during 
the  process  of  loading. 

The  main  macroscopic  effect  of  such  cracks  on  laminate  properties  is  reduction  of  stiffness  Since 
a  laminate  is  a  thin  plane  structural  element,  the  effect  of  interlaminar  cracks  on  stiffness  reduction  is 
in  general  negligible.  Indeed  the  only  effect  of  stiffness  reduction  of  such  cracks  appears  in  the  case  of 
bending  stiffness  at  regions  of  high  transverse  shear.  Here  we  shall  only  be  concerned  with  in-plane 
stiffness  and  in  this  case  only  the  intralaminar  cracks  produce  stiffness  reduction.  Thus  the  problem  to 
be  considered  is  the  effect  of  distributions  of  intralaminar  cracks  on  laminate  in-plane  stiffness. 

There  has  been  much  interest  in  this  problem  in  recent  years.  We  cite  the  works  of  Reifsnider[  1] 
and  Highsmith  and  Reifsnider[2],  who  treated  the  problem  in  terms  of  shear  lag  analysis.  Laws  and 
Dvorak[3],  who  employed  the  self-consistent  scheme  approximation  to  assess  the  stiffness  reduction 
of  a  cracked  lamina  in  a  laminate,  and  Talreja[4],  who  employed  continuous  damage  theory. 

The  various  approaches  cited  well-illustrate  the  dilemmas  with  which  one  is  faced  in  the  general 
problem  of  analysis  of  the  effects  of  damage  on  stiffness  reduction.  It  is  certainly  desirable  to  recognize 
and  incorporate  the  microstructural  details  of  damage  (cracks  in  the  present  case)  but  this  leads  to 
severe  mathematical  difficulties,  primarily  because  of  the  zero  traction  conditions  which  have  to  be 
satisfied  on  crack  surfaces.  Such  treatment  must  of  necessity  be  approximate,  thus  the  simple  shear 
lag  analysis.  It  is  thus  tempting  to  describe  damage  in  terms  of  some  continuous  fields,  as  was  done  in 
[4],  but  the  precise  mathematical  nature  of  such  a  description  is  not  clear,  must  be  assumed,  and 
involves  parameters  which  must  be  experimentally  determined  with  the  hope  that  they  are  indeed 
material  constants.  Here  we  present  a  microstructural  approach  based  on  variational  methods  which 
is  more  realistic  and  accurate  than  the  simple  shear  lag  analysis.  The  analysis  involves  construction  of 
admissible  stress  fields  which  satisfy  equilibrium,  all  boundary  conditions  on  laminate  and  crack 
surfaces  and  all  interlaminar  traction  continuity  conditions.  These  admissible  fields  deviate  from  the 
exact  stress  fields  only  in  terms  of  a  single  approximation  as  will  be  explained  below.  The  stress 
energies  associated  with  these  fields  define  lower  bounds  on  stiffness  which  are  optimized  by 
variational  methods. 

2.  ANALYSIS  OF  CRACKED  CROSS-PLY  UNDER  TENSION 

We  first  consider  a  cracked  cross-ply  laminate  which  is  subjected  to  simple  tension  (Fig.  1).  In 
the  absence  of  cracks  the  stresses  in  the  plies  are  plane  and  constant  in  each  lamina  and  are  readily 
found  by  elementary  conventional  laminate  analysis.  We  denote  such  stresses  as  a°yml  and 
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Fig.  1.  Cracked  cross-ply  loaded  by  tension  and  shear, 
where  m=  1,2  is  the  ply  index.  In  particular  we  shall  use  the  notation 


t<xi>  , 


*i. 


.0(2)  . 


where  from  now  on  1  indicates  the  90°  ply  and  2  indicates  the  0°  plies. 

The  presence  of  cracks  will  produce  a  different  space  variable  stress  system  which  is  xz  plane.  We 
denote  the  stress  perturbations  produced  by  the  cracks  in  the  different  plies  by  tr}"';  f,  j  =  x,  z.  We 
introduce  the  basic  approximation  that  a and  <x£>  are  not  functions  of  z.  We  thus  denote  these 
stresses  as 

<t"  - -<r, <£,(*)>  (2.1) 

where  </>,  and  <f>2  are  unknown  functions.  We  seek  stress  systems  e\}}  and  rr}^  which  satisfy  the 
equilibrium  equations  and  all  boundary  and  inteiface  conditions  of  the  laminate.  It  may  be  shown 
(see  [S]  for  details)  that  such  stress  systems  are 

-ff^x) 

<*,'  -  °\ <t>’(x)z 

-  (T.^fxfid.r,  -z1)  (2.2) 
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with  the  crack  boundary  conditions 

*<±«)«1.  *'(±«>-0  (2.3) 

where  t,  and  t2  are  ply  thicknesses,  <p  =  <t>{  and  primes  on  <t>  denote  x  derivatives.  The  stresses  are 
defined  for  any  region  of  length  la  between  adjacent  cracks  (Fig.  1)  and  these  distances  may  be 
different.  The  admissible  stresses  dtj  of  the  laminate  are  then 

+  7’.  (2.4) 
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The  stresses  (2.4)  satisfy  equilibrium,  traction  continuity  at  ply  interfaces  z  =  ±f|t  vanishing  of 
tractions  on  laminate  free  surfaces  z  =  ±h  and  on  crack  surfaces  x  *=  ±a  and  the  loading  conditions 
on  the  laminate  in  the  x  direction. 

Let  the  average  stress  in  the  x  direction  on  the  laminate  be  denoted  by  a0  thus, 

<x0  -  NxJ2h.  (2.5) 


It  is  then  rigorously  true  that  the  stress  energy  stored  in  the  laminate  is  given  by 


(2.6) 


where  £x  is  the  Young’s  modulus  in  the  x  direction  of  the  cracked  laminate.  The  complementary 
energy  functional  0C  of  the  cracked  laminate  is  given  by 

2  (2.7) 

where  SIJU  are  the  local  compliances  of  the  plies.  It  has  been  shown[5]  that  in  the  present  case  this  can 
be  expressed  as 


V+\  l  £  WJ'd''  (2-8) 

where  rr'"' are  given  by  (2.2)  and  £?  is  the  Young's  modulus  of  the  undamaged  laminate.  Introducing 
(2.2)  into  (2.8)  and  performing  z  integrations  we  have 


0 « 


<7q 
2  E°x 


(2.9) 


L. 
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j,  =  [coo^  +  Co^^  +  C^^y  +  Cn^yjdi  (2.10) 

where 

i  =  x/tu  p„  =  j./r, 

Coo  “  1/Ct"!" 

C02  -  (vt/£t)(a  +  2/3)-vax/3£a 
Cjj  =  (/l+  1)(3A2  +  12A  +  8)/60£T 
C„  =i<l/GT  +  A/GA) 

A  *  tj/t, 

where  a,  are  the  distances  between  adjacent  cracks  and  £A,  £T,  vA,  vT,  GT  and  GA  are,  respectively, 
axial  and  transverse  Young’s  moduli,  Poisson's  ratios  and  shear  moduli. 

It  follows  from  the  principle  of  minimum  complementary  energy  that 


or  in  terms  of  (2.6)  and  (2.9) 


l/,<  0t 


±<l+aj!iyj 


This  provides  a  lower  bound  on  Ex  for  any  4>n  which  satisfy  the  conditions  (2.3)  on  the  cracks.  To 
obtain  the  best  lower  bound  the  right  side  of  (2.12)  is  minimized  using  techniques  of  the  calculus  of 
variations.  This  procedure  establishes  differential  equations  for  the  functions  <t>H  which  all  have  the 
form 

dV,  dV 

d^ +pd?r+^--0 


p  —  (C02  —  CU)C22,  <?  —  C00/C2J 


with  the  boundary  conditions 


<M±p,)=  I-  ^-(±p»)  =  o 


P,  =  ajt ,. 


The  solution  of  (2.13)  subject  to  (2.14)  is 


<j>Ji)  =  A(ml)Cha£  cos  Pi  +  A'^Shai  sin  Pi 

a  =  q11*  cos  (6/2),  p  =  q11*  sin  (0/2)  (2.15) 

tan  0  =  v/(4q/p2  -  1 

where  -4J,1’  and  /tj,21  are  easily  determined  from  the  boundary  conditions  (2.14).  The  solution  is 
different  when  4 q/p1  <  1. 

Introducing  these  functions  into  (2.10)  and  (2.12)  we  find 

a'6’ 


UmiTWF*  VHTl mi ir»  W"w  v*  xrm  u-'m 
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where 

*i  -  ffi/ffo  ■  Vff° 

ff(A)  m  (3A2  +  l2A+8)/60,  A-tj/t,  (2.17) 

,  .  *-»/  t  ■  a*  Ch2ap,-cos Ifip, 

X(p.)  -  0(8  +  ft  )  ^  s.n  + psh2ap„ 

and  the  brackets  < — >  denote  average  with  respect  to  the  random  variable  p„.  In  the  case  when  the 
distances  between  cracks  are  equal  and  given  by  la 

<p>  -  p  -  a/r. 

F  1  (2.18) 
<z(p)>  -  z(p)- 

When  the  cracks  are  far  apart  (all  p„  are  large)  and  do  not  interact,  the  results  reduce  to 

Y^Jo+-£rkifl(XW*1+P1)tic  (2.19) 

where  c  is  the  crack  density,  i.e.  the  number  of  cracks  per  unit  length. 

When  distances  between  cracks  are  small,  thus  all  p,  ->0,  we  have 


—  <i-  _* L(±  J_\ 

Et"  E°x  +  1+\{et+ XEj 


Determination  of  the  Poisson's  ratio  vx>  (tension  in  the  x  direction,  strain  in  the  y  direction) 
needs  special  additional  analysis.  We  shall  consider  this  problem  in  approximate  fashion.  Tension  in 
the  y  direction  (Fig.  1)  will  primarily  produce  stress  a„  in  the  plies  and  very  small  stresses  axx.  If  the 
latter  are  neglected  it  follows  that  the  cracks  in  the  90°  ply  have  no  effect  on  the  stress  distribution 
and  therefore  E,  and  vyx  are  not  affected  by  the  cracks.  Thus 

Er-E° 

(2.21) 

V  -  V 

It  follows  from  the  usual  elastic  symmetry  of  compliances  that 


Therefore 


Em  E, 


V  ^  -  V 

"  £? V 


—  v° 
£0  v*r 


Experiments  to  determine  stiffness  reduction  of  glass/epoxy  cross-ply  laminates  have  been 
performed  by  Highsmith  and  Reifsnider  and  have  been  reported  in  [2].  The  elastic  properties  of  the 
single  unidirectionally  reinforced  ply  are 

£A  *=  41.7 GPa,  Et  =  13.0 GPa,  GA  =  3.40GPa 

CT  -  4.58  GPa,  vA  =  0.30,  vT  =  0.42 

where  GT  and  vT  have  been  computed  and  the  others  are  experimental  values.  For  the  [0  90}] 
laminate  it  follows  from  elementary  laminate  analysis  that 

E°  -  20.30  GPa 

k,  m  ojgo  =  0.636  (90°  ply). 
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Figure  2  shows  the  comparison  of  the  result  (2. 16M2. 18)  with  the  experimental  data  for  a  [0790 j], 
laminate  and  it  is  seen  that  the  agreement  is  excellent.  The  agreement  is  not  quite  as  good  for  the 
[0790°],  laminate  but  it  should  be  noted  that  in  this  case  there  is  an  unusual  discrepancy  between 
measured  (22.7  GPa)  and  calculated  ££  (27.57  GPa)  which  needs  explanation. 


X  ANALYSIS  OF  CRACKED  CROSS-PLY  UNDER  SHEAR 

The  cracked  cross-ply  (Fig.  1)  is  now  assumed  to  be  loaded  by  a  shear  membrane  force  Nir  This 
is  equivalent  to  average  shear  stress  loading  aMr  given  by 

“  Nx,l2h  «  t0.  (3.1) 

When  there  are  no  cracks  in  the  laminate  the  only  surviving  internal  stress  is 

<  -  (3.2) 

in  both  plies. 

When  the  90°  ply  is  cracked  as  in  Fig.  1,  the  state  of  stress  in  the  plies  is  antiplane,  thus  the  only 
surviving  stresses  are  a „  and  <x„.  We  proceed  as  in  Section  2  and  construct  an  admissible  stress 
system  on  the  basis  of  the  assumption  that  aXJ  is  a  function  of  x  only.  Thus 


'  ^-toCl-^M] 

»U'”To[l-^W]. 

The  functions  ^,(x)  and  i^2(x)  are  related  in  terms  of  the  equilibrium  condition 

*  To>»- 


(3.3) 


(3.4) 


Satisfaction  of  stress  equilibrium,  interface  traction  continuity  and  boundary  conditions  leads  to 
the  admissible  stresses 

*  t0[1-i^(x)] 


-  T<#(X)Z 
“  *07- *'(*)<*“-)■ 


(3.5) 


Proceeding  as  in  Section  2,  in  terms  of  the  principle  of  minimum  complementary  energy  we 
arrive  at  the  result 


"  l  +  <77yrp>/A//<p> 


:  _  3(1+  I/a) 

^  1+aGa/Gt’ 

Again  when  the  cracks  are  equidistant,  2a  apart. 


(3.6) 


*  =  t„t 


2/lf 


<P>  =  P  *  a/f , 
< Thup >  =  Thpp. 


(3.7) 
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YOUNG'S  ANO  SHEAR  MODULUS  REDUCTION 
OP  C0*/90\]t  GLASS/EPOXY  LAMINATE 


Fig.  3.  Young's  and  shear  modulus'reduction  of  [0790j],  glass/epoxy  laminate. 


When  the  distances  between  cracks  are  large 


x,>  1  +  1  /Amp 


and  when  the  distances  are  small 


Ga 

1+1/A 


(3.8) 


(3.9) 


The  last  result  has  a  very  simple  interpretation.  It  is  the  shear  modulus  of  the  cross-ply  when  the  shear 
stiffness  of  the  90°  ply  has  reduced  to  zero.  Figure  3  shows  the  reduction  of  shear  modulus  due  to 
crack  distribution  as  a  function  of  crack  density.  The  reduction  as  a  function  of  Young's  modulus 
given  in  Fig.  2  is  also  shown  by  comparison. 

4.  DISCUSSION  AND  CONCLUSIONS 

The  variational  method  developed  permits  the  evaluation  of  stiffness  reduction  of  any  [0m/90-], 
laminate.  Of  particular  interest  is  the  limiting  value  of  stiffness  when  the  cracks  become  very  close  to 
one  another.  Table  1  shows  such  results  for  glass/epoxy  and  graphite  (T300)/epoxy  quoting  the  ratio 
E,(p-*0)/E°x. 

It  is  seen  that  because  of  the  much  larger  ratio  EJEf  for  the  graphite/epoxy  the  stiffness 
reduction  is  much  smaller  since  the  undamaged  0°  plies  are  the  major  contribution  to  stiffness.  The 
situation  is  different  for  shear  modulus  where  the  limiting  ratio  is  defined  in  all  cases  by  the  simple 
expression  (3.9)  which  is  purely  geometrical.  The  analysis  also  easily  permits  evaluation  of 
approximate  stresses.  These  have  been  given  in  detail  in  [5]  and  their  significance  for  the  damage  to 
failure  process  of  the  laminate  has  also  been  discussed. 


Table  I 


Glass/epoxy 

Graphile/epoxy 

[0/90  1, 

0.770 

0.970 

[0-/90;], 

0.522 

0.914 

VR»i»/l*l 
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When  the  position  of  the  90°  and  0°  plies  is  reversed,  i.e.  the  0°  is  in  between  the  90°,  the  cracked 
laminate  is  of  course  easily  analyzed  by  the  same  method  and  it  may  be  shown  that  the  results  are  the 
same  as  for  the  90°  inside  and  the  0°  outside  provided  that  total  thickness  of  plies  of  each  kind  is 
preserved.  It  should  be  noted  that  the  laminates  analyzed  are  of  limited  practical  significance.  It 
would  be  desirable  to  evaluate  stiffness  reduction  in  cracked  [  ±  45°]„  [0°/  ±  45°],  and  [0°/90°/  ±  45°], 
laminates.  For  such  cases  a  minimal  requirement  is  to  analyze  the  [  ±  45°],  laminate  which  is  cracked 
in  both  layers,  thus  forming  an  orthogonal  crack  net.  This  is  of  course  much  more  difficult  than  the 
present  case  and  such  work  is  now  in  progress. 
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Analysis  of  Orthogonally  Cracked 
Laminates  Under  Tension 

The  problems  of  stiffness  reduction  and  stress  analysis  of  cross-ply  fiber  composite 
laminates,  where  all  plies  are  cracked  in  fiber  directions,  are  treated  by  a  variational 
method  on  the  basis  of  the  principle  of  minimum  complementary  energy.  The 
Young's  modulus  obtained  is  a  strict  lower  bound  but  is  expected  to  be  close  to  the 
true  value  on  the  basis  of  experience  with  a  previous  analysis.  Approximate  values  of 
Poisson 's  ratio  and  internal  stresses  have  been  obtained.  The  latter  reveal  important 
tendencies  of  continued  failure  by  delamination. 


Introduction 

The  major  damage  which  develops  in  laminates  under  static 
or  cyclic  loading  is  in  the  form  of  interlaminar  and  in¬ 
tralaminar  cracks.  The  former  develop  gradually  and  slowly  in 
between  plies.  The  latter  appear  suddenly  and  in  large 
numbers  in  plies  in  which  the  stresses  reach  critical  values, 
perhaps  defined  by  the  first  failure  criteria  of  the  plies.  They 
are  families  of  parallel  cracks  in  fiber  direction  and  their 
macroscopic  effect  is  reduction  of  the  inplane  stiffness  of  the 
laminate. 

The  subject  of  the  analysis  of  stiffness  reduction  of  cracked 
laminates  has  received  repeated  attention  but  has  mainly  been 
concerned  with  cross-ply  laminates,  i.c.,  [0^/90*],  configura¬ 
tions  in  which  only  the  90*  plies  are  cracked.  The  main 
methods  of  analysis  employed  are:  a  simple  shear-lag  method 
(Reifsnidcr  and  Talug,  1980;  Reifsnider  and  Jamison,  1982), 
self-consistent  approximation  to  assess  ply  stiffness  reduction 
in  conjunction  with  classical  laminate  analysis  (Laws  el  ai., 
1983,  1983),  and  a  variational  method  (Hashin,  1983).  All  of 
these  methods  give  results  which  are  in  good  to  excellent  agree¬ 
ment  with  experimental  data.  The  last  method,  unlike  the 
others,  also  provides  useful  estimates  of  the  internal  stresses  in 
the  cracked  laminate. 

In  many  cases  of  damage  in  laminates,  several  plies  will  be 
cracked  and  in  particular  adjacent  plies.  As  a  typical  and 
relatively  simple  case  we  shall  consider  the  problem  of  a 
(0,1/90*],  laminate  in  which  all  plies  have  intralaminar  cracks 
and  which  is  thus  orthogonally  cracked.  Because  of  the  com¬ 
plicated  interaction  of  the  orthogonal  cracks  this  problem  is 
much  more  difficult  than  the  one  with  plies  cracked  in  only 
one  direction.  We  are  aware  of  only  one  publication  in  the 
literature  (Highsmith  and  Reifsnider,  1986)  on  this  subject, 
which  is  concerned  with  evaluation  of  stresses  in  a  typical 
repeating  cracked  laminate  element  by  numerical  analysis. 
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Fig.  1  Orthogonally  cracked  laminate 


The  present  problem  is  closely  related  to  that  of  a  cracked 
±43*  symmetric  laminate  which  is  indeed  merely  a  rotated 
0*/90*.  It  is  easily  seen  that  uniaxial  tension  of  a 
l+45j,/-45*],  in  bisector  direction  can  be  analyzed  in  terms 
equibiaxial  tension  and  shear  of  a  (0* /90*],  laminate.  Indeed 
a  typical  case  of  orthogonal  cracking  is  encountered  for  cyclic 
loading  of  a  ±43*  laminate.  The  problem  of  shearing  of  a 
cross-ply  will  be  considered  elsewhere.  Here  we  are  concerned 
with  uniaxial  tensile  loading  of  an  orthogonally  cracked  cross- 
ply.  Our  purpose  is  to  evaluate  stiffness  reduction  due  to 
cracks  and  approximate  local  stresses.  We  shall  do  this  by 
generalization  of  the  variational  method  which  has  been 
developed  in  Hashin  (1983). 


872 1  Vol.  54,  DECEMBER  1987 


Transactions  of  the  ASME 


[MUUUUURT 


my  -gtnrJir  ^"WWgWTOt1; iru i 


i»i4»i4ii4»*- 

— h  ■("  h  — ■ 


U CRACKED 


Fig.  2  (Updating  •tamanl 

Admissible  Field  Construction 

Consider  an  orthogonally  cracked  laminate,  Fig.  I,  and  a 
typical  laminate  element  defined  by  intersecting  crack  pairs, 
Fig.  2.  Let  the  laminate  be  subjected  to  constant  membrane 
force 

(1) 

and  no  other  loads.  Thus  a*  is  the  average  stress  6„  over  the 
laminate  and  also  over  the  thickness  2A. 

In  the  event  that  longitudinal  and  transverse  crack  families 
are  each  equidistant,  the  stresses  and  strains  in  each  element 
(Fig.  2)  are  the  same,  the  boundary  and  interface  conditions  of 
the  repeating  element  can  be  formulated  and  the  problem  can 
be  treated  numerically  in  terms  of  Finite  elements,  if  desired. 
The  problems  associated  with  such  an  approach  will  be 
discussed  at  the  end  of  the  paper.  Here  we  shall  proceed  dif¬ 
ferently;  it  is  our  purpose  to  construct  admissible  stress  Fields 
which  satisfy  equilibrium  and  all  traction,  boundary  and  inter¬ 
face  conditions.  These  admissible  Fields  will  then  be  optimized 
in  the  context  of  the  principle  of  minimum  complementary 
energy  to  yield  approximate  stresses  and  lower  bounds  on 
stiffness. 

The  stresses  in  the  undamaged  laminate  under  the  loading 
considered  are 

=  o*“  =  a<;>=*<V 

o«»-o®-*SJ>  o*  o*»  -o®«  *<?>«•  l2) 

where  the  k  coefficients  are  easily  found  from  conventional 
laminate  analysis. 

Let  9„  be  a  three-dimensional  admissible  stress  Field  within 
the  cracked  laminate.  By  definition  such  a  stress  field  must 
satisfy  equilibrium,  traction  continuity,  and  traction  bound¬ 
ary  conditions.  For  reasons  of  symmetry  it  is  sufficient  to  con¬ 
sider  one  half  of  the  element  of  Fig.  2  defined  by  -  a  .  x  » 
a,  -b  m  y  *  b,  0  m  Z  =  h,  as  will  be  understood  from  now 
on.  We  list  required  boundary  and  interface  conditions 

3iiW,0}-0  a”W.0)-0  (3) 

*0  *0  (4) 

agW.M-0 

(5) 
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aJ!?(  ±ajj)  -  a™(  ±aj,i)  -  a*i>(  ±aj\z)  *  o 

(6) 

3<«(x,  ± bj)  -  a®  (x,  ±b,z)  -  d®(x.  ± bj)  -  0 

Here  and  from  now  on  superscripts  1,  2  indicate  the  plies, 
equations  (3)  are  symmetry  conditions,  equations  (4)  are  free 
surface  conditions,  equations  (3)  express  traction  interface 
continuity,  and  equations  (6)  are  zero  traction  conditions  on 
the  crack  surface.  There  are  still  needed  traction  conditions  on 
the  faces  x  »  ± a,  /,  £  z  =  h  and  y  «  ±b,  0  £  z  £  f, 
which  shall  be  considered  further  below. 

The  xy  in-plane  parts  of  the  admissible  stress  Fields  are 
chosen  in  the  form 

atf-oi'UI -♦,(*>!  a®-c?>[l-*j(x))  (7a) 

d£W,l,(l-^,0>)l  a®=o®[i-*20’)l  (7  b) 

(7c) 

where  and  are  unknown  functions.  The  physical 
significance  of  these  assumptions  will  be  discussed  further 
below. 

Force  equilibrium  of  the  undamaged  laminate  in  the  x  and  y 
directions  requires 

(8) 

«<;»/,+ e®/, -o 

Examining  force  equilibrium  in  the  same  directions  in  the 
cracked  laminate,  in  terms  of  the  admissible  stresses  (7),  we 
find 

«i"f|*i(*)  +  dPh+iW  *  0 

This  reduces  the  number  of  unknown  functions  from  four  to 
two. 

In  view  of  equation  (7c),  the  equilibrium  equations  for  the 
admissible  ply  stresses  assume  the  form 

9UJ,  +  *«.*“  0 

+  =  0  (10) 

where  commas  denote  partial  differentiation. 

We  now  derine  perturbation  stresses  o'"1  for  the  plies,  m 
■  1.2,  by 


+oi"'  (ID 

where  oj'*'  are  given  by  the  undamaged  laminate  stresses  (2). 
Insertion  of  equations  (7)  into  (10)  for  each  ply  and  systematic 
integration,  with  elimination  of  unknown  residual  functions 
by  use  of  equations  (3)-(3),  and  utilization  of  equations  (9), 
leads  to  the  results 


off-  -oi'»«(x) 

(12a) 

*2— —  oi*»P(x) 

02/) 

(*)z 

02  b) 

eg1  -  '  (x)  ( h — z ) 

02  g) 

’)« 

(12 d) 

<4  O’)  (*-«) 
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*U>-  W'**  U) +<>^Ck)1-j-(A/1  -z2)  (12c) 

•ff— J-M'9'  (x)  +  <>  foU-fA-z)*  (IV) 

where  prime  and  dot  indicate  x  and  >  derivative,  respectively, 
and 


♦  **i 

X*/j/r| 

It  follows  from  equations  (6)  and  (9)  that  d  and  must 
satisfy  the  boundary  conditions 

4(±a)«l  ♦'(*<»)«  0 

*(±a)«l  t(±b)-0 

The  above  formulation  determines  the  traction  components 
on  the  faces x  -  ±a,tl  .  z  «  band>>  »  ±b,  0  ^  Z  ^  ti¬ 
lt  is  seen  that  the  only  nonvanishing  traction  component  is 


(13) 


(14) 


a2(±fl)*oi'»(l+/,/tj) 


US) 


Thus  the  admissible  tractions  are  continuous  across  all  com¬ 
mon  faces  of  adjacent  elements,  regardless  of  their  a,  b  dimen¬ 
sions.  Therefore,  the  stresses  (11)-(12)  rigorously  satisfy  all 
conditions  of  admissibility. 


Variational  Formulation 

Next  we  evaluate  the  complementary  energy  functional  0C 
associated  with  the  admissible  stresses.  In  the  present  case 
where  tractions  are  prescribed  everywhere 


U'—^S^vdV 


(16) 


where  V  is  laminate  volume  and  S„w  are  local  compliances.  It 
has  been  shown  in  Hashin  (1983)  that  for  any  cracked  body 


where 


utm’T 


O'c—Y^oro-udV 


•o5  +  «ri 


U,$0C 

Combining  equations  (17)  and  (19)— (21)  we  have 


1 


1  21/; 

-  + - T- 


(22) 


(17) 

(18o) 
(180) 
(18c) 

Here  a?  are  the  actual  stresses  in  the  uncracked  body,  i.e.,  in 
the  undamaged  laminate;  thus  U ?  is  the  actual  stress  energy  of 
the  undamaged  laminate.  In  the  present  case  it  is  rigorously 
true  that 

--ar’’  <'» 

where  £J  is  the  Young's  modulus  in  the  x  direction  of  the  un¬ 
damaged  laminate.  Furthermore,  the  stress  energy  Ue  of  the 
cracked  laminate  can  be  rigorously  expressed  in  the  form 
a*1 

U'—^r-V  (20) 

where  £,  is  the  effective  Young’s  modulus  of  the  cracked 
laminate.  Also  from  the  principle  of  minimum  complementary 
energy 


E,  £?  Va°‘ 
which  provides  a  lower  bound  on  Ex. 

In  the  present  case  the  stresses  (12)  are  identified  with  the  a'v 
in  equation  (18c).  In  order  to  evaluate  equation  (18b)  we  list 
the  local  stress  energy  densities  in  the  plies  as  referred  to  the 
common  x,  y,  z  system  of  the  laminate 

2W»-d"f/EA  +  (o»>1  +o<i»1)/£r— 

+a^)2rA/EA -<t"<t"2*T/ET 

+  «>*  +  VCA  +  otf  /CT  (23) 

2#**  -  dg/EA  +  («?/  +c2>J)/£r-a«>(o«' 

+  o™)2*a  /Ea  -  a^y^2eT(ET 

♦(•gf+offVOa+ogf/Gr 

where  the  coefficients  are  in  terms  of  the  following  properties 
of  the  unidirectional  ply  material 

Ea  *  Axial  Young's  modulus  (fiber  direction) 

¥a  =  Axial  Poisson’s  ratio 
£r*  Transverse  Young’s  modulus 
rr*  Transverse  Poisson’s  ratio 
CA  »  Axial  shear  modulus 
Cr*  Transverse  shear  modulus 

Suppose  a  typical  element.  Fig.  2,  has  in-plane  dimensions 
2 am  and  2 b„.  Then  for  half  the  laminate  0  =  z  =  b 

Oi-p  f"  P  Wdxdydz 

+  f  (  f  HA2>dxdydz  (24 a) 

0'C  =  Z  U-„  (24 b) 

which  can  now  be  evaluated  in  terms  of  equations  (12).  For 
simplicity  in  writing  we  shall  perform  detailed  analysis  for  the 
case  when  am  »  a,  b„  »  b  and  give  results  for  the  more 
general  case  of  unequal  intercrack  distances  at  the  end.  When 
am  and  b„  are  constant  we  can,  without  loss  of  generality,  in¬ 
terpret  all  preceding  energies  as  those  stored  in  an  element 
with  dimensions  2a,  2b,  h.  Introducing  equations  (12)  into 
(24a),  performing  all  z  integrations  and  using  the  nondimcn- 
sional  variables 

( =x/t,  t|  =>>//,  (25) 

we  have 

2 0;-«‘>,/}P'  P  [klA^  +  lk,k,BM 

j -#i  J  ~n 

+  k2C9+2+Alk2*'t  +  Blk2t2 

+A1kIb(k,<k‘  +  k,+)+B1k,Hk,<t>’  +k,i) 

+  C(M*+*^)Jld(d, 

where  prime  and  dot  now  denote  (  and  i\  derivative,  respec¬ 
tively,  and 


(26) 


(21) 


♦  -♦(() 
Pi  ■«/»! 

k,  -*<•' 


i  ='Hv) 

Pi  =  b/t, 
k, 


(27) 
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(28) 


At  *  I /\EA  +  UEt  B0  *  -  (1  +  1/\)>a/Ea 

c0*i/ea  +  i/\et 

A ,  *(k/CA  +  l/Gr)/3  B,  *(l/GA  +  \/Gr)/3 

a2  -  l(3X + 2)pt/Et  -  \pA  /£a  )/3 
BJ  **  K3X  +  2)» a  'Ea  ~  X»r/£rJ/3 

C-  (X+  IK3X1+  12X  +  8)/60£r 
and  the  boundary  conditions  (14)  now  assume  the  form 
*(±p,)«‘l  *'(±p,)«  0 

^(±p,)-l  ^(±p2)*0 

In  view  of  (22)  any  functions  4  and  ^  which  satisfy  equa¬ 
tions  (29)  will,  when  introduced  into  equations  (24)  and  (22), 
provide  an  upper  bound  for  1  !EX.  To  obtain  the  best  bound 
U'c  should  be  minimized.  Using  techniques  of  the  calculus  of 
variations  we  find  that  the  minimum  conditions  of  the  integral 
(26)  are  the  integro-differeniial  equations 

d*+ 


(29) 


dV 

rfV 

where 


+Pi 


+P» 


k. 

B0  1  f-i 

+9i*+ 

k. 

C  2p2J-^ 

k. 

B0  1  [»* 

dr,*  +  **  + 

k. 

C  2p,J-„ 

A}  — A 1 

\ 

Aq 

*"  C 

<' ~C 

B2-B, 

Co 

c 

*-c 

-rj— (  <t>dt=4>  y— ( 

2pi  J  -#i  2pj  J  -#j 


m. 


0=<M£) -"*i^ 
*  =  * o(’i)-tn2<fi 
.  Bo 


**  B0 


Xj,  Aq 


a  *  ql/4cos(0/2) 


tan  g»V4q/pJ-  I 
when4q<pJ  p<0 
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f(v)»Chav  g(i>)  *  Chfiv 


(36) 

V(  —  1  ±Vl  -4q/pi)p/2 

where  v  is  either  (  or  tj,  and  p  and  q  are  either  of  p, ,  p2  or  q, , 
9s- 

The  case  p>0  is  unlikely  for  the  usual  stiff  fiber  composites. 

To  solve  the  problem  (29)  and  (30)  let  equations  (34)  be  ex¬ 
pressed  in  the  form 

*(l)=.D,/i  ({)+*•, 8,  ({)-/"■* 

Ml)  =£>1/1(1)  +F1g1(i\) -m2* 
where/, ,  g,  indicate  solutions  in  terms  of  p, ,  q,  and/2 ,  g2 — in 
terms  of  p2,  g2.  Introducing  equations  (37)  into  (29)  we  obtain 
sets  of  linear  equations  for  the  constants  D,  F  with  solutions 

(l+ffii^)g,(Pi) 


Z),= 


F  i  =  ~ 


/i<Pi)8i(Pi)-/i(Pi)8i(Pi) 

(l+m,^) 


(30o) 

(30b) 

(31) 


/i(Pi)«iO>i)'/i(Pi)8i(Pi) 
(l+/n2*)g2(p2) 


(38) 


/s  (Pi)8i  (Ps)  (Ps)8s(Ps) 


F,=  - 


(1  +m2j)/,(p,) 
/2(P2)«s(Ps)  -?s(Ps)8s(Ps) 


Now  average  both  sides  of  equations  (37)  as  in  (33).  The  result 
can  be  written  in  the  form 


where 


Solution  of  the  Integra-Differential  Equations 

Obviously  the  solution  of  equations  (30)  with  (29)  is  insen¬ 
sitive  to  change  of  sign  of  independent  variables.  Therefore, 

*(£)=*<-«)  *(lW(~1>  (32) 

Now  define  the  mean  values 

r  n 


“r 


^  +  +Ftgi  *  (1  +/n, $)u, 

£  +  m2d  *=  A/2  +  ^282  *  (1  +/n20)u2 
/.8,'(P,)-/,'(P,)8, 


(40) 


/i(Pi)8i'(Pi)-/i'(Pi)8i(Pi) 


«2  = 


/i(Pj)82(P2) 


(41) 


(33) 


Since  expressions  (41)  are  known,  equations  (40)  can  be  solved 
for  i  and  Thus 


Then  the  general  solutions  of  equations  (30)  can  be  written  in 
the  form 


u,  —  m ,  ( 1  -  w,)w2 
1  -m,/n2(l  -u,)(l  -u2) 


(34a) 

(34b) 

(34c) 


j  u}-m2(l-u3)<J, 

1  —  m,m2(l  “  w,)(l  -  wj) 


(42) 


where  the  terms  with  subscripts  zero  are  general  solutions  of 
the  homogeneous  versions  of  equations  (30)  and  the  remaining 
terms  are  particular  solutions.  A  homogeneous  differential 
equation  of  type  (30)  has  four  independent  solutions  which 
may  be  arranged  into  two  symmetric  and  two  antisymmetric 
ones.  In  view  of  equations  (32)  only  symmetric  solutions  need 
be  retained.  These  are: 

when  4q>p1 

/(h)*  ChavcosQv  g(v)*  Shavsin&v 
0  *  ql/4sin(0/2) 


These  results  together  with  equations  (30)-(39)  determine  the 
solution  (37)  of  equations  (29)  and  (30).  Specific  results  for  (he 
two  cases  (33  and  (36)  are  given  in  the  Appendix. 

The  perturbation  stresses  (12)  now  assume  the  form 


(35) 


<$}=  -*<»*({) 


o JJI»  -«i0(Mi|) 


eg!-—  c<r"\Mi) 


(43) 


$ 


I 

t, 
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Effective  Young's  Modulus  £,  and  Poisson’s  Ratio  wx. 

It  is  first  necessary  to  evaluate  the  minimum  value  of  the  in¬ 
tegral  (26)  and  to  introduce  it  into  (22).  Direct  evaluation  of 
the  integral  (26)  in  terms  of  the  functions  4  and  i  which  have 
been  determined  above  is  a  formidable  undertaking  which  can 
fortunately  be  avoided.  The  best  way  to  proceed  is  to  integrate 
equation  (30a)  multiplied  by  4  from  -p,  to  p,  and,  similarly, 
equation  (30b)  multiplied  by  i  from  -pj  to  p2.  These  in¬ 
tegrals  obviously  vanish  and  when  transformed  by  integration 
by  part  with  use  of  the  boundary  conditions  (29)  they  can  be 
identified  with  substantial  parts  of  the  integral  (26).  As  a  result 
of  this  procedure  we  find  the  simple  result 

(p,)  +  JkJp,iftpi)]  (44) 

To  evaluate  the  third  derivatives  in  equation  (44)  we 
evaluate  the  mean  of  equation  (30a)  with  respect  to  (  within 
l— Pi.Pil  and  of  equation  (30b)  with  respect  to  i)  within  [  -p2, 
Pj].  Noting  that  the  means  of  the  second  derivatives  vanish 
because  of  equations  (29)  and  taking  into  account  equations 
(31)  we  obtain 

-C±l^-A0*+-£-fl0* 

P| 

(45) 

Pj 

Substituting  equations  (45)  into  (44)  and  (27)  we  have 

U'Cmm  m  2/1abo°2  KV10  +  kfBa)k,4>  +  (kxB0  +  k,C0)k,t]  (46) 

Thus  the  minimum  energy  has  been  expressed  in  terms  of  $ 
and  i  which  are  explicitly  determined  by  equations  (42).  Since 
the  reference  volume  in  the  present  case  is 

K*4abb-4ab(/,+rj)  (47) 

we  have  from  equations  (46),  (47),  and  (22) 

X  ~  £f +TTx(^+/r^) 


*  (kjA#  +  k,Bt)kx  Kj  m  (kxBa  +  krCa)kf 

This  provides  a  rigorous  lower  bound  on  E,  in  the  event  that 
all  cracks  are  equidistant.  If  this  is  not  the  case  let  inter-crack 
distances  be  2 am  in  x  direction  and  2b,  in  y  direction,  respec¬ 
tively.  Then  the  admissible  stress  fields  employed  are  still  valid 
in  each  laminate  element  of  dimensions  2 am,  2 b„  h.  The  com¬ 
plementary  energy  functional  is  now  equations  (24)  and  each 
term  in  the  summation  is  expressed  in  terms  of  functions 
and  in  the  form  equation  (26),  leading  to  extremum  condi¬ 
tions  of  form  (30)  for  each  of  these  functions,  subject  to 
boundary  condition  of  type  (29) 

*..(  ±P|»)  M  1  *i(±Pi«)=0 


This  defines  means  and  which  are  precisely  4  and  i  ex¬ 
pressed  in  terms  of  (30).  Thus 

qi* <51> 

and  the  relevant  laminate  volume  is  now 

Y-4b£  £a„b,  (52) 

Introducing  equations  (St)  and  (52)  into  (22)  we  obtain  again  a 
lower  bound  for  E, : 

This  result  can  be  put  into  more  compact  form  if  the  inter¬ 
crack  distances  are  random  variables  a,  b,  assuming  values  am , 
b„.  Then  (50)  are  the  values  of  associated  random  variables  p, 
and  p2  with  joint  probability  density  fuoction  P(plt  p2).  In 
this  event  (22)  will  assume  the  form 

1  <  1  Kx<0lp14(pl)>  +K,<piPi\l/(pi)> 

E,  m  Ei  (1+X)<p,p2> 

where  <  >  indicates  probability  average.  It  can  probably  be 

assumed  that  p,  and  p}  are  independent  random  variables.  In 
this  event 

J,(Pi.b:)  =  /,,0>.)'>2<P2)  (54) 

and  (S3)  assumes  the  form 


_L  <  ±+-±-\K 

Ex  ~  E*  1  +  xl** 


<«Mp,)>  „  <\Mp2)> 


«» 


Pi»  Pj ,-b,//, 


<Pi>*j  <Pj>  (56) 

and  also 

<p,>*<o>/r,  <p2>  *  <b>//2 

The  effective  Poisson’s  ratio  »x,  is  defined  by 

4r-  (57) 

'if 

Since  vXf  cannot  be  bounded  by  variational  methods  we 
evaluate  it  in  terms  of  the  approximate  stresses  resulting  from 
our  variational  treatment.  For  the  case  of  equidistant  cracks  in 
the  x  and  y  direction 

f->  !•  !,*  (")dxdydz} 


It  follows  from  the  ply  orientations  and  the  stress-strain  rela¬ 
tions  of  the  ply  material  that 

‘iV  “  “  /Ea  +  <*»/ea  ~  va  /Ea 


*®=  “  ^a  +  ag>/£  r  -  o£'  »t/Et 

Next  we  introduce  the  stresses  (11)  and  (12)  into  equations  (59) 
and  the  resulting  expressions  into  equation  (58).  We  observe 
that:  (I)  The  stresses  a°(m|  produce  the  Poisson’s  ratio  v°x/  of 
the  undamaged  laminate;  (2)  The  contribution  of  au  to  equa¬ 
tion  (58)  vanishes  since 

T-.  J*^'(jlf)‘ted>»2bj%'{x)(ix  =  2b[«'(x)l*.. 

and  this  vanishes  because  of  equations  (14).  A  similar  result  is 
true  for  the  integral  of  +  over  the  rectangle.  An  easy  calcula¬ 
tion  then  yields  the  result 
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*  Table  1  Properties  of  unidirectional  material.  Moduli  in 
GPa. 

_ ea  et  Ga  Gt  K a  "t 

Class/Epoxy  41.7  13.0  3.40  4.38  0.300  0.420 

♦  Graphite/ Epoxy  208.3  6.5  1.63  2.30  0.255  0.413 


o  Jk-  (  E*  E *  \ 

**  £«  “  1+X  V  £r  “  Ea  ) 


(60) 


where  i  is  defined  by  equations  (33).  Evaluation  of  vyx  is  of 
course  entirely  analogous. 


Results,  Discussion  and  Conclusions 

The  approximate  nature  of  the  analysis  which  has  been 
developed  here  is  contained  in  the  assumed  form  of  the 
stresses  (7)  which  certainly  cannot  be  correct.  Indeed  from  ex¬ 
perience  with  crack  field  solutions  one  would  expect  local 
singularities  at  the  crack  tips,  which  are  absent  from  equations 
(7).  But  it  should  be  remembered  that  in  the  usual  structural 
laminates  we  are  dealing  with  very  thin  plies  of  thickness  of 
order  0.2-0.4  mm.  Typical  glass  and  carbon  fiber  diameters 
are  0.01  mm.  Assuming  for  the  sake  of  argument  that  fibers 
are  arranged  in  regular  hexagonal  arrays,  we  find  that  at  the 
usual  0.60  fiber  volume  fraction  the  ply  thickness  can  accom¬ 
modate  about  16-32  fibers.  Now  the  basic  assumption  in 
laminate  theory  is  that  plies  are  homogeneous  anisotropic  with 
elastic  properties  equal  to  the  effective  elastic  properties  of  the 
unidirectionally  reinforced  ply  material.  But  it  should  be 
remembered  that  the  concept  of  effective  modulus  is  based  on 
stress  and  strain  averages  over  representative  volume  elements 
(RVE)  which  by  definition  must  contain  many  fibers.  It 
follows  that  a  ply  can  accommodate  perhaps  2-3  RVE  sizes 
across  its  thickness.  Therefore,  a  severe  gradient  of  local 
average  stress  or  strain  through  the  thickness  is  impossible. 
This  does  not  present  a  problem  for  the  simple  cases  of  un¬ 
damaged  laminates  when  the  stresses  are  constant  or  linear 
across  ply  thickness,  but  the  situation  is  very  different  in  the 
case  of  an  intralaminar  crack.  If  it  is  assumed  that  the  fields  of 
the  homogeneous  ply  model  are  averages  over  RVE,  then  the 
crack  tip  singularities  of  fracture  mechanics  and  their 
associated  large  stress  gradient  cannot  exist.  If  this  assumption 
is  not  made,  there  remains  the  choice  between  two  difficult 
alternatives:  (a)  The  crack  tip  is  in  the  matrix  and  recognizes 
the  adjacent  fibers  as  distinct  heterogeneities,  thus  creating  an 
intractable  problem  in  crack  mechanics;  (b)  The  crack  tip 
vicinity  is  governed  by  a  nonlocal  elasticity  theory  of  unknown 
nature  or  with  many  unknown  elastic  constants.  None  of  these 
alternatives  is  useful  or  practical  and  it  would  therefore  appear 
that  the  simple  assumptions  (7),  with  their  resultant  gentle 
through-the-thickness  variations  of  the  stresses,  are  less  severe 
than  would  perhaps  appear  at  first  sight. 

All  of  these  considerations  apply  also  to  finite  element 
analysis.  The  number  of  finite  elements  allowable  through  the 
thickness  must  be  severely  limited  and  singular  elements  raise 
the  same  problems  as  the  crack  tip  singularities  discussed 
above. 

Highsmith  and  Reifsnidcr  (1986)  have  treated  the  problem 
of  orthogonally  cracked  cross-plies  in  terms  of  a  refined 
laminate  theory  advanced  by  Pagano  (1978)  in  which  stress 
variation  through  ply  thickness  is  restricted  to  be  linear.  The 
treatment  is  of  necessity  restricted  to  equidistant  cracks  and 
requires  extensive  numerical  analysis.  Their  results,  Highsmith 
and  Reifsnider  (1986)  show  that  the  stress  oa  depends  only 
weakly  on  y,  the  stress  a„  depends  only  weakly  on  x,  and  a„  is 
negligible.  This  is  in  accordance  with  our  basic  assumptions 
(7).  Furthermore,  the  shear  stress  a„  practically  does  not  de¬ 
pend  on  y  and  aa  practically  does  not  depend  on  x  and  this 
agrees  with  (12  b.dj.i). 
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For  purpose  of  numerical  evaluation  we  have  used 
glass/epoxy  and  graphite/epoxy  laminates  with  unidirectional 
ply  properties  as  in  Hashin  (1985).  These  are  given  in  Table  I . 

Table  2  shows  Young’s  modulus  reduction  and  also 
Poisson’s  ratio  for  various  laminates  with  equidistant  cracks 
in  both  directions.  Also  shown  by  comparison  are  results  of 
Hashin  (1985)  for  the  case  when  there  are  cracks  only  in  the 
90’  ply.  It  is  seen  that  the  effect  of  the  cracks  in  the  0*  plies  is 
quite  small,  no  doubt  because  of  the  small  value  of  transverse 
a„  stress  in  the  0*  plies  due  to  tension  in  0*.  i.e.,  x  direction. 
This  is  also  in  agreement  with  the  results  of  Highsmith  and 
Reifsnider  (1986).  The  situation  is  quite  different  for  the 
Poisson’s  ratio.  It  is  seen  that  cracks  have  significant  effects 
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Fig.  3  #2|,  «2I>  *»4  #21  variation 


on  the  value  of  this  quantity  and  that  the  effect  of  orthogonal 
cracks  can  be  very  substantial  as  is  evidenced  in  particular  by 
the  Graphite/Epoxy  cross  ply. 

For  the  loading  case  of  uniaxial  tension  considered  here  the 
stresses  of  major  interest  are:  and  eg?  which  are  transverse 

to  the  Fibers  in  their  respective  layers  and  may  thus  produce  in¬ 
tralaminar  cracks:  the  interlaminar  shear  stresses  a„  and  o„', 
the  stress  <su  normal  to  the  plane  of  the  laminate  which  is  call¬ 
ed  peeling  stresses  when  tensile.  It  is  seen  from  the  form  of 
equations  (12)  and  (43)  that  the  shear  stresses  assume  their 
maximum  values  at  the  interface  of  the  plies  and  that  oa 
assumes  its  maximum  value  on  the  midplane  z  ■  0.  However , 
the  value  of  au  at  the  interface  is  also  of  importance,  for  the 
stress  which  produces  interlaminar  separation  may  be  smaller 
than  the  one  producing  midplane  transverse  cracking. 

We  have  evaluated  stress  distributipns  for  (0V90*), 
graphite/epoxy  laminates  assuming  equidistant  cracks  a  *  b, 
for  various  crack  densities  defined  by  different  ratios  p  = 
a//, .  We  have  observed  the  following  significant  phenomena: 

(1)  The  tensile  stress  a?)  is  much  smaller  than  the  tensile 
stress  •“  the  load  direction.  It  can,  therefore,  be  assumed 
that  during  loading  no  new  intralaminar  cracks  will  develop  in 
the  0*  plies  in  the  y  direction. 

(2)  The  stresses  o2?  and  °a  w  close  in  value  to  the  same 
stresses  produced  in  the  same  laminate  when  only  the  90  plies 
are  cracked  (Hashin,  1983).  A  comparison  is  shown  in  Fig.  3 
for  the  case  «/r»  -  2. 

(3)  The  maximum  value  of  tensile  transverse  stress  o2/  is 
midway  between  cracks,  is  equal  to  the  corresponding  stress 

in  the  undamaged  laminate  when  distances  between  cracks 
are  large,  and  decreases  monotonically  with  increasing  crack 
density  (decreasing  p)  (Fig.  4). 

(4)  The  intralaminar  shear  stress  an  has  significant  values 

(Fig.  3)  and  therefore  the  important  interlaminar  shear  stress 
is  r  ■  .  Figure  4  shows  the  value  of  as  a  func 

lion  of  crack  density.  It  is  seen  that  this  stress  starts  out  with  a 

878/ Vol.  54,  DECEMBER  1987 


value  of  0.48a2>  which  is  certaintly  significant.  When  the 
vlaue  of  is  large  enough  to  cause  cracking,  r  may  be  large 
enough  to  produce  shear  delamination.  It  is  interesting  to  note 
that  increasing  crack  density  decreases  this  shear  stress. 

(3)  The  peeling  stress  aa  assumes  significant  values. 
Figure  3  shows  the  variation  of  midplane  oa  over  one  quarter 
of  the  square  defined  by  the  intersection  of  two  pairs  of 
cracks.  Observe  the  large  tensile  stress  1 .08a2*  in  the  middle  of 
they  *  a  edge.  When  oi,)  is  large  enough  to  crack  the  90*  ply 
transversely,  this  stress  is  certain  to  produce  a  crack  in  the  xy 
plane.  Note  that  for  the  present  laminate  au  at  the  ply  inter¬ 
face  is  one  half  of  the  midplane  value.  Depending  on  interface 
strength  such  a  stress  may  also  produce  interface  separation. 
The  largest  value  of  au  is  always  located  on  the  edges  y  =  ±a. 
For  small  crack  density  the  maximum  is  close  to  the  cracks  x 
=  ±a  and  its  location  moves  to  the  center  of  the  edge  with  in¬ 
creasing  crack  density.  Figure  4  shows  the  variation  of  max  o„ 
with  intercrack  distance.  It  would  be  of  interest  to  analyze 
various  other  cases  and  also  biaxial  loading,  which  is  easily 
carried  out  by  superposing  the  results  of  two  uniaxial  loads  in 
the  x  and  y  directions. 
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appendix 

Denote  for  convenience 

^  =  +j 

For  case  (35) 

A  -2(1  +/»,*, )[a,C/t(a,pl)sintf,p1) 
+  /JtSA(a,P|)cos(/31pl)]/Af| 

F"i  -2(1  +m,*j)(/3,CA(a,p,)sin03|P,) 
-a,5A(a,p,)cosOSlpl)J/A/l 
Mk  -a,sin(2^,pl)  +  0l5A(2a,p,) 

u  -  ^CACfratpQ-cos^p,)! 

‘  p,(a,//J,+0,/«i)Af| 


To  obtain  Dlt  Fit  and  ui2,  simply  replace  I  by  2  in  the  above 
equations.  Then 

♦({)“D,CA(a,f)cos03l{)+F,SA(ol()sin03lO-m,^ 

Ml)  «DjCA(aji|)cos03j,)  +F,SA(a2n)sin(/32n)  -m20 
For  case  (36) 

D,  -  (I  +  m.+jW.SAOJ.p.VN, 

f|  -  “(I +m,+2)a,SA(o,p,)/N, 

Af,  -^|CA(a|P,)5A(^,pl)-a|SA(olP|)CA(^lpl) 

«i  -Oi/«i  -or,//S,)SA(or1pl)SA(/J,pl)/plA/l 

To  obtain  Dlt  Flt  and  u2,  replace  1  by  2  in  the  above  equa¬ 
tions.  Then 

+({)«0,Gft(a,{)+^|SA(0,{)-ml{ 

^(l )  *D2CA(or2ij)  +FlShtfJii)  -m20 
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Abstract.-  The  various  damage  phenomena  in  statically  or  cyclically  loaded 
unidirectional  fiber  composites  and  their  laminates  are  described  and 
classified.  Existing  analytical  methods  to  evaluate  the  effect  of  damage  in 
laminates  on  their  thermoelastic  properties  are  described  ar.d  discussed. 
Fundamental  problems  arising  in  application  of  classical  mechanics  methods  to 
very  thin  cracked  plies  are  pointed  out.  Existing  work  on  damage  prediction 
and  failure  analysis  is  described. 

1.  Introduction 

The  word  Damage  has  different  connotations  in  different  areas  of  science 
and  technology,  not  to  speak  of  everyday  life.  In  the  present  context  we 
shall  define  damage  as  an  accumulation  of  many  small  internal  defects  inside  a 
solid  body.  Defect  implies  an  internal  discontinuity  such  as  void,  crack  or 
dislocation.  it  is  assumed  that  the  damaged  body  can  be  described  on  a 
suitable  scale  of  magnitude  by  an  effective  continuum  with  stress-strain 
relations  defined  a3  the  relations  between  average  stress  and  average  strain 
where  averages  are  taken  over  the  usual  repr esentati ve  volume  elements  (RVE) 
which  are  small  compared  to  body  dimensions  yet  contain  a  large  number  of 


defects; 
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The  mechanics  of  a  damaged  body  will  be  called  Damage  Mechanics.  We 
shall  distinguish  between  micro-damage  mechanics  (MIDM)  and  macro-damage 
mechanics  (MADM).  The  former  implies  evaluation  of  the  detailed  micro-fields 
of  stress  and  strain  in  the  vicinity  of  defects  which  may  then  be  utilized  to 
evaluate  average  stress  and  strain,  thus  the  effective  or  macro  stress-strain 
relations,  and  also  to  obtain  some  information  concerning  local  failures. 
Within  this  frame  the  defects  are  explicitly  recognized  by  shape  and  spatial 

V. 

distribution.  By  contrast  the  latter  approach,  often  also  called  continuum 
theory  of  damage,  describes  damage  in  terms  of  abstract  fields  of  scalar, 
vectorial  or  tensorial  nature  which  also  contain  unknown  parameters  to  be 
determined  by  fit  to  experimental  data.  Such  description  of  damage  is  by  no 
means  unique  and  the  appropriate  choice  may  at  this  time  be  regarded  as  an 
open  question.  Evidently  MADM  cannot  be  expected  to  provide  information 
beyond  the  dependence  of  macro-material  parameters  or  effective  stress-strain 
relations  of  the  damaged  body  on  the  chosen  "damage"  field. 

It  would  appear  therefore  that  MIDM  should  be  preferred  since  it  is  more 
closely  tied  to  physical  reality  and  because  it  provides  more  relevant 
information.  It  is  unfortunately  also  very  difficult  to  carry  out  since  it 
requires  the  determination  of  detailed  fields  produced  by  many  interacting 
defects  and  this  can  be  done  only  for  relatively  simple  geometries  or  on  the 
basis  of  simplifying  assumptions.  It  is  therefore  our  view  that  MIDM  and  MADM 
should  be  viewed  as  complementary  approaches.  Whenever  MIDM  is  possible  it 
3hould  be  pursued;  otherwise  we  should  resort  to  MADM.  It  is  of  primary 
importance  that  MIDM  and  MADM  should  communicate  efficiently  and  share  their 
results.  A  recent  IITTAM  Symposium  on  Mechanics  of  Damage  and  Fatigue  [3]  has 
been  concerned  with  damage  mechanics  from  both  those  points  of  view. 


Damage  In  the  composite  materials  which  are  used  in  engineering  practice 
occurs  primarily  in  the  form  of  micro-cracks.  There  exists  a  considerable 
body  of  literature  on  stiffness  reduction  of  homogeneous  elastic  bodies  by 
crack  distributions  most  of  which  is  based  on  what  has  been  here  called  MIDM 
and  well  illustrates  the  difficulty  of  the  problems  involved.  Such 
difficulties  are  magnified  in  composites  because  of  the  presence  of  a 
reinforcing  phase  in  the  form  of  inclusions  or  fibers  and  consequently  cracks 
do  not  only  interact  with  one  another  but  also  with  the  lnhomogeneities.v 

We  shall  here  focus  our  attention  on  unidirectional  fiber  composites  (UFC) 
and  laminates  which  are  made  of  UFC  plies  with  analytical  emphasis  of  the 
latter.  From  the  engineering  point  of  view  there  are  two  main  reasons  to 
study  damage  in  fiber  composite  laminates.  Firstly  the  state  of  damage  is  a 
prelude  to  the  event  of  failure.  Several  attempts  to  predict  laminate  failure 
in  simple  fashion,  without  consideration  of  damage,  can  be  found  in  the  older 
literature  and  these  undertakings  have  not  been  very  successful.  It  would 
therefore  appear  that  analysis  of  damage  is  prerequisite  to  analysis  of 
failure. 

Secondly,  it  is  necessary  to  have  an  estimate  of  the  remaining  fatigue 
life  and  residual  strength  of  a  damaged  laminate  and  to  characterize  the  state 
of  damage  by  observable  parameters.  In  metal  structures  this  is  done  in  terms 
of  the  present  and  (predicted)  future  of  a  dominant  crack.  But  in  laminates 
damage  develops  as  a  collection  of  many  interior  microcracks  and  therefore 
different  methods  are  needed.  It  might  be  hoped  that  stiffness  reduction, 
measured  and  analyzed,  could  serve  as  such  parameters;  hence  the  fundamental 
importance  of  analytical  relation  of  stiffness  to  crack  distribution  or 
density. 
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2.  Classification  of  Damage 

A  UFC  is  made  or  stiff  and  strong  parallel  continuous  fibers  (glass, 
graphite,  carbon)  which  are  embedded  in  a  relatively  compliant  matrix  (polymer, 
aluminum).  We  are  here  primarily  concerned  with  the  case  of  polymeric  matrix. 
When  such  a  UFC  is  loaded  by  tension  in  fiber  direction,  whether  static  or 
cyclic,  fibers  will  rupture  at  different  random  locations,  Fig.  1.  These  fiber 
ruptures  are  micro-cracks  which  may  be  regarded  as  damage  in  the  sense  defined 
in  the  introduction  above.  At  each  rupture  there  must  develop  a  significant 
shear  interfacial  stress  as  may  be  easily  shown  ju3t  from  equilibrium 
considerations,  Rosen  [27].  This  shear  stress  may  produce  interfacial 
deterioration  or  separation  in  the  form  of  interfacial  cracks  which  is  another 
kind  of,  subsequent,  damage.  A  statistical  analysis  of  strength  based  on  this 
kind  of  damage  has  been  pioneered  in  [27]  and  has  been  refined  and  extended  in 
subsequent  works,  e.g.  [21]. 

IT  on  the  other  hand  the  usual  thin  specimen  of  UFC  is  loaded 
transversely,  normal  to  the  fibers,  the  specimen  fails  abruptly  by  a  single 
crack  along  the  fibers  without  prior  discernible  development  of  internal 
damage.  There  are  thu3  two  major  modes  of  failure  in  tension:  The  fiber  mode 
which  Involves  prior  development  of  damage  and  the  matrix  mode  which  does  not, 
[7]. 

It  is  of  considerable  Interest  to  evaluate  the  stiffness  reduction  of  a 
UFC  which  has  been  damaged  by  tension  in  fiber  direction.  Steif  [28]  has  given 
an  approximate  analysis  of  reduction  of  Young's  modulus  E^  in  fiber  direction, 
on  the  basis  of  a  simple  shear  lag  model.  Reduction  of  transverse  modulus  E^ 
and  axial  shear  modulus  by  cracks  along  fibers  has  been  analyzed, 

approximately,  on  the  basis  of  so-called  self-consistent  estimates  in  [6]  and 
[17].  All  of  these  treatments  fall  within  the  MIDM  category. 


5 


Since  the  strength  and  stiffness  of  UFC  are  very  high  in  fiber  direction 
but  low  transversely  to  the  fibers  it  i3  necessary  in  most  engineering 
applications  to  use  these  materials  in  the  form  of  laminates  as  is  shown 
schematically  in  Fig.  2.  We  consider  for  the  sake  of  simplicity  only  laminates 
which  are  symmetric  with  respect  to  their  mtdplanes  and  are  subjected  to 
constant  midplane  membrane  loads.  Then,  as  is  well  known,  the  state  of  stress 
in  any  ply  of  the  laminate  is  constant  and  plane  except  for  edge  boundary 

V 

layers  where  the  state  of  stress  is  three  dimensional  and  complex. 

Static  or  cyclic  loading  may  lead  to  the  following  events: 

(a)  If  there  is  high  stress  in  fiber  direction  a  ply  will  develop  fiber  damage 
as  has  been  described  above.  This  will  be  called:  intralaminar  fiber  damage. 

(b)  At  some  level  of  transverse  tensile  stress  a  fiber  mode  crack  will  be 
triggered  by  some  micro-defect.  Such  an  intralaminar  crack  propagates  along 
the  fibers  and  generally  reaches  the  edge3  of  the  laminate.  The  formation  of 
such  a  crack  is  a  random  event  and  the  probability  of  its  occurrence  increases 
with  increasing  magnitude  of  tensile  transverse  stress.  It  follows  that, 
Initially,  only  a  small  number  of  isolated  intralaminar  cracks  will  appear  in 
the  ply.  Isolated  Implies  that  there  is  no  mutual  crack  interaction  and 
therefore  the  transverse  tensile  stress  builds  up  from  zero  on  the  crack  to 
its  initial  magnitude  when  there  are  no  cracks.  This  is  unavoidably 
accompanied  by  Interlaminar  shear  stress,  Fig.  3,  a  phenomenon  similar  to  that 
of  interfacial  shear  stress  near  a  fiber  rupture.  As  the  ply  tensile  stress 
increases  because  of  increase  of  external  load  more  and  more  intralaminar 
cracks  will  appear.  Fig.  M  redrawn  from  X-radiographs  given  in  [30]  shows  the 
increase  of  crack  density  in  the  90°  ply  of  a  graph!  te/ epoxy  cross-ply 
laminate,  with  load  increase. 
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As  cracks  become  closer  the  In-between  tensile  stresses  diminish  and  can 
no  longer  build  up  as  before  and  therefore  ever  increasing  rate  of  loading  is 
required  to  produce  new  cracks.  This  is  well  illustrated  by  plots  of  crack 
density  versus  load  as  shown  in  [133*  Experimental  evidence  shows  that  crack 
densities  in  plies  reach  asymptotic  saturation  values  at  which  they  form 
roughly  regular  patterns  [2*1,25].  This  has  been  called  Characteristic  Damage 
State  (CDS),  Reif snider  [2*1]. 

\ 

We  shall  call  the  type  of  damage  described  Intralaminar  matrix  damage. 

(c)  With  intralaminar  saturation  the  interlaminar  stresses  produced  near  the 
cracks,  and  In  particular  at  locations  where  intralaminar  cracks  in  neighboring 
plies  cross,  may  produce  local  interfacial  debonding  between  plies  [26]  which 
may  be  described  as  a  collection  of  small  interlaminar  cracks  or  disbonds. 
This  will  be  called  interlaminar  damage. 

Note  that  a  laminate  may  have  Initial  interlaminar  damage,  before  loading, 
in  the  form  of  interlaminar  disbonds,  due,  for  example,  to  faulty 
manufacturing.  This  damage,  by  itself,  ha3  no  effect  in  the  case  of  in-plane 
loading  of  a  laminate  and  only  little  effect  in  the  case  of  bending.  The 
reason  is  that  in  the  former  case  the  plies  are  in  states  of  plane  stress; 
therefore  the  ply  interfaces  are  traction  free  and  thus  an  interlaminar 
disbond,  sufficiently  far  removed  from  the  edges,  has  no  effect.  In  the  latter 
case  there  are  also  Interlaminar  shear  stresses  due  to  bending  which  are,  in 
principle,  affected  by  such  a  disbond  but  are  in  general  very  small  compared 
to  ln-plane  stresses. 

In  the  case  of  ln-plane  compressive  load  a  disbond  whose  dimensions  are 
much  larger  than  the  smaller  distance  from  its  plane  to  laminate  free  surface 
may  be  the  source  of  local  buckling  and  propagation  of  an  interlaminar  crack. 
This  is  an  isolated  damage  phenomenon  and  is  therefore  not  within  our  present 
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(d)  It  has  been  observed  [15,16]  that  intralaminar  cracks  trigger  the 
development  of  fiber  ruptures  in  adjacent  plies,  which  are  distributed  in  bands 
along  the  cracks. 

(e)  Independently  of  the  damage  modes  described,  the  high  interlaminar  stress 
concentration  at  laminate  free  edge  may  be  the  source  of  initiation  and  growth 

of  an  interlaminar  edge  crack  which  may  propagate  into  the  laminate  with 

\ 

increasing  load  or  number  of  load  cycles.  The  analysis  of  such  a  crack  is 
within  the  scope  of  fracture  mechanics  and  since  It  is  an  Isolated  damage 
event,  such  as  the  disbond  buckling  mentioned  above,  it  is  not  within  the  scope 
of  our  discussion  of  damage. 

Most  of  the  analytical  discussion  to  follow  will  be  concerned  with 
intralaminar  matrix  damage.  It  is  seen  from  our  description  of  damage  that  it 
is  the  central  damage  mode,  as  it  is  the  source  of  others.  It  is  the  primary 
cause  of  stiffness  reduction  of  a  laminate  and  most  of  the  analytical 
investigations  of  laminate  damage  have  been  concerned  with  this  damage  mode. 

3.  Damage  Mechanics  of  Laminates 

3-1.  Stiffness  reduction  and  stress  analysis 

According  to  the  description  and  discussion  of  the  various  forms  of 
damage  in  laminates  given  above  the  primary  problem  is  to  evaluate  the  effect 
of  distributions  of  many  Intralaminar  cracks  in  the  plies.  The  major  problems 
are:  (a)  Analysis  of  stiffness  reduction,  (b)  Analysis  of  internal  stresses  and 
consequent  Implications  for  the  Internal  failure  process. 

All  of  laminate  analysis  in  the  literature  is  based  on  the  assumption 
that  the  pile's  are  homogeneous  anisotropic  with  the  effective  properties  of 
the  UFC.  This  is  a  valid  assumption  for  undamaged  laminates  for  the  usual 
cases  of  membrane  force  and  moment  loading  when  the  variation  of  the  primary, 


in-plane,  laminate  stresses  through  ply  thickness  is  constant  or  linear.  But 
the  presence  of  intralaminar  cracks  may  invalidate  this  assumption  because  it 
would  lead  to  very  large  stress  gradients  through  very  thin  plies.  The 

Implications  of  this  difficulty  will  be  discussed  further  below.  Even  the 
classical  version  of  the  problem  is  exceedingly  difficult  for  it  requires  three 
dimensional  analysis  of  the  effects  of  many  interacting  cracks  within  different 
plies  inside  a  laminate. 

We  shall  first  summarize  various  simple  approaches  which  have'  been 
developed.  Within  the  context  of  an  old  and  somewhat  primitive  attempt  to 
estimate  laminate  failure  it  has  been  assumed  that  ply  failure  is  governed  by 
some  semi- empirical  failure  criterion,  e.g.  a  quadratic  polynomial  of  stress, 
thus  determining  first  ply  failure  which  mostly  occurs  in  the  matrix  mode. 
Since,  as  has  been  explained  above,  the  failed  ply  will  develop  many 

intralaminar  cracks,  it  is  assumed  that  the  transverse  Young's  modulus  E^.  and 
the  axial  shear  modulus  of  the  "failed"  ply  reduce  to  zero.  A  new  stress 
analysis  of  the  laminate  is  then  performed  on  the  basis  of  stiffness  reduction 
of  the  failed  ply  and  the  next  ply  to  fail  is  identified  and  its  stiffnesses 
are  then  similarly  reduced.  This  procedure  which  has  sometimes  been  called 
the  Ply  Discount  Method  Is  continued  until  the  occurrence  of  fiber  failure  in 
the  0°  ply,  i.e.  the  one  with  fiber  orientation  in  load  direction.  This  simple 
method  is  based  on  a  crude  estimate  of  the  effect  of  intralaminar  damage  and 

not  surprisingly  such  analysis  of  failure  is  mostly  unreliable. 

Reifsnider  [2*i]  and  Relfsnlder  et  al.  [25,26]  have  advanced  beyond  this 
simple  approach  by  recognizing  that  when  a  ply  develops  intralaminar  cracks 
the  load  transf er  must  take  place  through  interlaminar  shear  stresses  in  the 
vicinity  of  the  crack.  Fig.  3-  They  have  accordingly  devised  a  simple 
"str  ength-of-materi  al  s"  type  shear  lag  analysis  which  recognizes  only  two 
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kinds  of  stress:  load  carrying  ply  stresses  which  are  constant  through  ply 
thickness  and  variable  in  load  direction,  and  interlaminar  shear  stresses  which 
are  confined  to  a  narrow  layer  of  arbitrarily  assumed  thickness  in-between 
plies.  While  the  local  stresses  predicted  by  such  simple  analysis  can  hardly 
be  reliable  the  stiffness  reduction  thus  determined  is  often  from  fair  to  good 
agreement  with  experimental  data. 

Another  approach  has  been  taken  by  Laws  and  Dvorak  [5],  They  have  used 

•  \ 

the  approximate  method  known  as  the  Self  Consistent  Scheme  (SCS)  to  estimate 
the  stiffness  reduction  of  a  cracked  ply  [17]  and  then  proceeded  on  that  basis 
to  evaluate  stiffness  reduction  of  cracked  laminates.  The  SCS  has  first  been 
employed  to  evaluate  stiffness  reduction  of  a  homogeneous  body  containing  many 
cracks  by  Budlansky  and  O'Connell  [M],  It  is  based  on  the  drastic  assumption 
that  any  crack  is  embedded  in  the  effective  medium  whose  properties  are  to  be 
determined.  The  results  thus  obtained  are  not  always  acceptable.  Thus  for 
randomly  oriented  penny-shaped  cracks  the  stiffness  vanishes  at  some  definite 
crack  density  which  is  unreasonable  for  stiffness  must  vanish  asymptotically 
with  Indefinitely  increasing  crack  density. 

The  difficulties  are  compounded  in  the  case  of  a  cracked  ply.  Fig.  5a 
shows  the  Crack  Opening  Displacements  (COD)  of  intralaminar  cracks  which  are 
of  necessity  restrained  by  the  neighboring  plies.  The  SCS  cannot  cope  with 
this  situation.  It  need3  to  assume  that  ply  thickness  is  Infinitely  larger 
than  typical  crack  length  a3  is  shown  in  Fig.  5b.  It  follows  that  (a)  crack 
density  of  Fig.  5a  has  to  be  somehow  translated  into  a  crack  density  parameter 
appropriate  to  Fig.  5b.  (b)  Crack  lengths  are  now  of  the  order  of  fiber 
diameters  and  In  spite  of  this  it  is  assumed  that  cracks  are  embedded  In  a 
material  with  the  effective  properties  of  the  UFC. 
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With  all  of  these  drastic  assumptions  the  SCS  stiffness  predictions  are  in 
fair  to  good  agreement  with  some  available  experimental  data  [5  3.  Obviously 
the  SCS  cannot  provide  any  Information  regarding  microstresses  and  local 
fail  ure. 

The  methods  described  so  far  may  be  regarded  as  approximate  MIDM 
approaches.  A  MADM  approach  has  been  proposed  by  Talreja  [29].  He  has 

described  interlaminar  and  intralaminar  damage  in  terms  of  general  vector 

\ 

fields  and  has  reduced  them  to  more  specific  forms  on  the  basis  of  symmetry 
arguments.  The  final  forms  involve  undetermined  constants  to  be  determined  by 
fitting  to  experiment.  The  basic  question  which  arises  is:  are  the  constants 
so  determined  material  constants?  i.e.  independent  of  laminate  Internal 
geometry.  To  examine  this  the  prediction  of  stiffness  reduction  due  to 
intralaminar  damage  in  terms  of  such  continuum  damage  fields  was  compared  to 
experimental  data  thus  obtaining  a  set  of  constants  for  glass  epoxy.  These 
constants  were  then  used  to  compare  analytical  to  other  experimental  results, 
obtaining  fair  to  good  agreement  [29].  However,  the  analysis  of  interlaminar 
damage  is  unacceptable  for  it  predicts  In-plane  stiffness  reduction  due  to 
Interlaminar  cracks  alone  and  as  has  been  pointed  out  in  part  2  cracks  have  no 
effect  on  in-plane  stiffness  (this  is  an  exact  elasticity  theory  conclusion). 

It  should  be  realized  that  stiffness  reduction  analysis  is  much  less 
sensitive  to  simplifying  assumptions  than  stress  analysis,  Tor  effective 
stiffness  is  defined  as  a  relation  between  stress  and  strain  averages  and  an 
average  is  not  necessarily  sensitive  to  local  errors.  This  may  explain  why 
3uch  simple  methods  often  yield  effective  elastic  properties  which  are  in 
reasonable  agreement  with  experimental  results. 

More  accurate  methods  are  clearly  required.  Bearing  in  mind  that 


analytical  rigorous  solutions  can  not  be  constructed  there  remain  two  avenues 
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of  approach:  More  accurate  analytical  approximations  and  numerical  analysis 

in  terms  of  finite  elements  (or  finite  differences).  All  existing  work  of  this 
nature  is  based  on  the  classical  description  of  a  laminate  in  terms  of 
homogeneous  anisotropic  plies.  Before  proceeding  we  must  examine  the 
implications  of  the  homogeneous  ply  assumption  within  the  present  context. 

A  typical  ply  thickness  may  be  of  the  order  of  .2  mm  while  typical  fiber 

diameter  (glass,  carbon,  graphite)  is  of  the  order  of  .01  mm.  At  typical  fiber 

\ 

volume  fraction  .60  such  ply  thickness  can  accommodate  about  15  fibers.  Now 
the  concept  of  effective  media  is  based  on  scaling  involving  three  orders  of 
magnitude,  [8],  the  microscale,  the  mlniscale  and  the  macroscale.  In  a 
composite  material  micro  is  the  size  of  heterogeneities  (e.g.  fiber  diameter). 
Mini  (al so  called  Meso)  is  the  size  of  a  typical  Representative  Volume  Element 
(RVE)  which  is  the  smallest  element  of  the  composite  to  exhibit  the  effective 
properties  with  sufficient  accuracy  and  must  contain  many  heterogeneities 
(fibers).  Macro  i3  the  size  of  the  composite  material  body.  To  apply 
elasticity  theory  in  terms  of  effective  elastic  moduli  to  a  composite  material 
body  there  must  be  fulfilled  the  necessary  requirement  Micro  <<  Mini  <<  Macro. 
Note  that  in  the  classical  theory  of  elasticity  there  is  no  microscale:  the 
differential  element  is  the  miniscale  while  the  body  dimensions  are  the 
macroscale. 

Furthermore,  in  classical  elasticity  there  are  only  surface  forces  and  no 
surface  moments  on  the  faces  of  a  differential  element.  Such  assumption  can 
be  transferred  to  RVE  faces  in  terms  of  local  surface  averages  only  if 
gradients  of  local  average  stress  are  sufficiently  small,  otherwise  non- 
classical  elasticity  theory  -  such  as  multipolar  or  non-local  -  would  have  to 
be  invoked  for  the  continuum  description  of  the  composite. 
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All  of  this  has  significant  implications  for  crack  mechanics  of 
composites.  Consider  as  a  pertinent  example  a  transverse  crack  through  a  UFC. 
It  must  be  our  primary  goal  to  avoid  the  intractable  problem  of  a  crack  which 
"can  distinguish"  between  fibers  and  matrix  and  therefore  the  UFC  must  be 
viewed  as  a  homogeneous  material  with  the  effective  properties.  A  typical 
plane  RVE  for  random  fiber  locations  is  a  square  of  about  10x10  fibers.  The 
crack  length  is  the  macro  dimension  which  must  be  by  order  of  magnitude  larger 
than  the  RVE  dimension.  If  we  view  this  crack  in  the  classical  crack  mechanics 
sense  then  the  usual  stress  "at  a  point"  is  now  replaced  by  local  average 
stress  and  this  stress  variable  will  become  singular  at  the  crack  tip.  This 
really  implies  that  local  average  stress  at  this  location  has  some  large 
unknown  value  and  therefore  a  criterion  for  crack  criticality  has  been 
formulated  in  terms  of  surface  energy,  this  being  regarded  as  a  material 
property.  The  contribution  of  the  crack  to  the  energy  balance  is  evaluated  in 
terms  of  the  stress  field  or  the  COD  in  spite  of  the  crack  tip  state  of 
uncertainty;  obviously,  such  evaluation  can  be  reliable  only  when  the  crack  tip 
uncertainty  regions  are  small  relative  to  total  crack  length  and  not  when  they 
are  substantial  parts  of  crack  length. 

In  the  case  of  an  intralaminar  crack  the  crack  length  is  equal  to  the  ply 
thickness  which  is  of  the  order  of  a  bout  two  RVE  sizes.  Since  this  is  the 
macro  dimension  the  necessary  condition  Min  <<  Macro  is  thus  violated.  Or  to 
put  is  differently,  the  crack  consists  entirely  of  crack  tip  regions  and  is  too 
short  to  develop  COD  and  stress  which  can  be  reliably  predicted  by  classical 
crack  mechanics  in  terms  of  effective  elastic  UFC  properties. 

This  is  a  severe  problem  which  goes  to  the  fundaments  of  mechanics  of 
solids.  It  implies  that  very  accurate  numerical  solutions  involving  of 
necessity  many  finite  elements  through  the  thickness  are  not  necessarily 


correct  and  singular  finite  crack  tip  elements  are  not  necessarily  applicable 
for  it  is  unlikely  that  there  are  crack  tip  singularities  of  local  RVE  average 
stress.  There  is  of  course  the  microfield  alternative  to  recognize  fibers  and 
matrix  at  the  crack  tip  but  this  is  an  intractable  problem,  even  from  the 
point  of  view  of  numerical  analysis,  not  to  mention  its  futility. 

We  shall  refer  to  the  problem  discussed  above  as  the  thin  ply  problem  and 

we  believe  that  it  is  unresolved  at  the  present  time  and  should  be  pursued. 

\ 

In  the  meanwhile  there  is  no  choice  but  to  continue  with  damage  analysis  based 
on  classically  homogeneous  plies,  allowing  however  certain  flexibility  in 
result  Interpretation  in  view  of  the  thin  ply  problem. 

The  most  basic  problems  concern  a  cracked  laminate  which  is  subjected  to 
constant  membrane  forces  which  produce  the  average  stresses 

0^  -  N  /2h  id  -  1.2  (1) 
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Here  E,,  E2,  are  Young's  moduli,  v,2,  v2,  are  Poisson's  ratios,  and  C12  is  the 

\ 

shear  modulus  -  all  for  the  cracked  laminate.  For  the  uncracked  laminate 
(2-3)  retain  their  forms  and  we  shall  distinguish  the  properties  in  this  case 
by  superscript  0 e.g.  E°  . 

It  is  seen  from  (2)  that  to  determine  the  effective  thermoelastic 
properties  of  a  cracked  laminate  it  is  necessary  to  evaluate  the  average 
strains  for  specified  load.  This  requires  in  general  a  detailed  internal  field 
solution,  satisfying  zero  traction  conditions  on  crack  surfaces  and  other 
boundary  and  continuity  conditions. 

The  effective  elastic  compliances  may  also  be  expressed  in  terms  of  the 
crack  opening  displacements  (COD),  [8].  Thus  for  applied  0.. 
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where  the  surface  integral  is  again  over  all  cracks  and  T®  is  the 
traction  at  crack  surface  location  in  the  uncracked  laminate. 

All  the  relations  listed  so  far  are  exact  and  may  be  used  to  obtain  the 
compliances  of  a  cracked  laminate  on  the  basis  of  exact  analytical  or  accurate 
numerical  solutions. 

An  important  and  more  flexible  alternative  is  the  variational  formulation, 
primarily  in  terms  of  the  principle  of  minimum  complementary  energy,  which  has 
been  initiated  for  cracked  laminates  in  [9].  Let  be  an  admissible  stress 
field.  Such  stress  field  must  satisfy  (a)  Equilibrium  everywhere,  (b)  Traction 
continuity  at  ply  Interfaces.  (c)  Laminate  load  boundary  conditions  and  crack 
surface  zero  traction  conditions.  In  this  case  the  compl ementary  energy 
functional  0°  is  defined  as 
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All  of  the  variational  work  and,  to  the  best  of  our  knowledge,  all  of  the 
numerical  work  has  been  concerned  with  the  relatively  simple  case  of  cracked 
CO®,  90®]s  laminates,  also  called  cross-ply  laminates.  Initially  the  case  when 
only  the  90°  ply  is  cracked,  Fig.  6,  has  been  considered  and  more  recently  also 
the  case  when  all  plle3  are  cracked.  Fig.  7.  These  will  be  referred  to  as 
singly  and  orthogonally  cracked,  respectively.  The  first  case  occuhs  when  a 
cross-ply  is  loaded  monotonically  or  cyclically  in  simple  tension,  the  second  - 
when  a  which  is  merely  a  rotated  cross-ply,  is  loaded  in  simple 

tension  and  also  when  a  cross-ply  is  subjected  to  cyclic  temperature  change. 

Basic  to  all  analytical  work  is  the  Identification  of  a  typical  or 
repeating  element  of  the  cracked  laminate.  Such  an  element  Is  indicated  in 
Fig.  7  and  is  shown  enlarged  in  Fig.  8.  In  the  event  that  in  each  crack  family 
the  cracks  are  equidistant  the  element  shown  is  a  repeating  element  in  which 
case  stress  and  strain  fields  are  periodic  with  periods  2a  and  2b.  Such 
periodicity  is  essential  for  numerical  work.  If  the  cracks  are  not  equidistant 
the  element  is  called  typical  and  it  is  then  a  typical  building  block  of  the 
cracked  laminate. 

An  admissible  stress  field  for  an  orthogonally  cracked  laminate.  Fig.  7, 

under  tension,  has  been  constructed  in  [10]  starting  out  from  the  premise  that 

o  is  only  a  function  of  x,  o  is  only  a  function  of  y  and  a  «  0.  Then 
xx  yy  xy 

systematic  integration  of  the  equations  of  equilibrium,  taking  into  account 
traction  continuity  and  boundary  conditions  defines  an  admissible  stress  field 
in  terms  of  two  unknown  functions  <K x )t  y ) ,  as  follows: 
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where  superscripts  1  and  2  identify  the  plies  as  labeled  in  Figs.  7,  8, 


o(2),  and  o(2)  are  ply  stresses  in  the  uncracked  laminate,  prime  and 


dot  indicate  x  and  y  derivative,  respectively,  and  the  functions  <t>  and  must 


satisfy  the  following  boundary  conditions 
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Note  that  the  stress  fields  thus  defined  satisfy  traction  continuity  where  the 


typical  elements  are  Joined,  thus  satisfying  admissibility  for  the  entire 
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cracked  laminate.  It  is  seen  that  in-plane  stresses  are  constant,  shear 
stresses  linear  and  o ^  is  quadratic  through  ply  thickness,  thus  exhibiting 
none  of  the  classical  crack  tip  singularities.  In  view  of  the  discussion  of 
the  thin  ply  problem  given  above  such  gentle  through-ply  thickness  variation 
does  not  seem  unreasonable. 

To  employ  the  variational  formulation,  (12)  is  introduced  into  (7).  Note 
in  this  respect  that  the  plies  are  transversely  isotropic  unidirectional  fiber 
composites  and  therefore  the  integrand  in  (7)  for  a  ply  has  the  form  v 
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where  1  indicates  fiber  direction  and  2,3  are  transverse  directions.  These 
material  axes  change  orientation  from  ply  to  ply.  Subscript  A  indicates 
property  in  fiber  direction  and  T  in  direction  transverse  to  fibers.  Carrying 
out  the  z  integration  in  (7)  we  obtain  a  double  integral  in  x,y  which  must  be 
minimized  to  obtain  an  optimum  bound.  This  results  in  linear  integro- 
differential  equations  for  the  function  $  and  ij»  subject  to  the  boundary 
conditions  (13).  These  equations  have  been  solved  analytically  and  then  (9)  or 
(11)  define  optimum  lower  bounds  on  Young's  modulus.  Furthermore  the 
optimizing  function  $  and  i|>  now  determine  approximate  stresses.  Details  are 
given  in  [10].  Note  that  Poisson's  ratios  cannot  be  bounded  directly,  but  they 
may  be  evaluated  approximately  by  averaging  of  the  strains  associated  with  the 
approximate  stresses. 

The  case  when  only  the  90°  ply  is  cracked  is  of  course  a  special  case  of 
orthogonal  cracking  when  the  Intercrack  distances  b,  say,  become  very  large. 
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A  detailed  solution  has  been  given  in  [9]  which  also  includes  results  for  shear 
modulus.  Note  that  results  for  shear  modulus  of  orthogonally  cracked  cross- 
plies  are  not  available  in  the  literature. 

Fig.  9  shows  analytical  results  for  reduction  of  Young's  modulus  of 
singly  cracked  [0#/90#  ]  glass/epoxy,  [93,  and  experimental  results  for  this 

®  3 

case  from  [133,  demonstrating  excellent  agreement.  Fig.  10  shows  stresses 
near  an  isolated  crack  in  90°  ply  for  [0°/90°3a  T300/epoxy  and  Fig.  11  for  the 
same  material  when  the  cracks  are  close  enough  to  interact  significantly.  It 
is  seen  that  in*  the  first  case  the  tensile  stress  axx  builds  up  to  its  initial 
value  in  the  uncracked  laminate.  But  in  the  second  case  it  cannot  do  so 
because  of  crack  proximity.  Therefore  a  significant  portion  of  ply  tensile 
force  must  now  be  carried  by  interlaminar  shear  stress  and  this  will 
ultimately  result  in  delamination.  Also  note  the  tensile  interlaminar  stress 
at  mid-distance  between  cracks  which  may  also  produce  delamination. 

Young's  moduli  of  orthogonally  cracked  laminates  are  only  insignificantly 
affected  by  the  crack  family  in  load  direction.  Thu3  to  all  practical  purposes 
they  are  the  same  as  in  the  previous  case  when  only  the  90°  ply  is  cracked. 
The  Poisson's  ratios,  however,  are  significantly  different  in  the  two  cases  of 
one  or  two  crack  families  as  is  shown  in  Fig.  12.  Fig.  13  shows  maximum 
values  of  in-plane  tensile  stress  <jxx,  principal  interlaminar  shear  stress  and 
of  tensile  o  on  [0V900]  laminate  midplane  as  function  of  crack  density.  It 

ZZ  3 

is  assumed  that  all  cracks  are  equidistant  and  that  crack  density  is  the  same 

in  both  plies.  The  ply  material  is  T300/epoxy  and  stress  values  are 

normalized  by  the  tensile  o  ply  stress  of  the  uncracked  laminate.  Note  the 

xx 

large  value  attained  by  a  which  is  bound  to  produce  cracking  on  the  laminate 
midplane.  Such  cracks  are  similar  to  Interlaminar  damage  discussed  in  Section 


2. 
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Analysis  of  singly  cracked  cross-plies  [9]  has  shown  that  elastic 
properties  approach  asymptotic  limits  with  increasing  crack  density.  Such 
limits  indicated  by  superscripts  S,  are  shown  below 


Material  Glass/Epoxy  Graphite/Epoxy 


Configuration 

[0/90 ]g 

[°/90j]3 

[0/90 ]g 

[0/90,]s 

E*7E° 
x  xy 

.770 

.522 

.970 

.91  4 

G*  /° 
xy  xy 

.500 

.250 

.500 

.250 

Note  that  Young's  modulus  reduction  of  graphite/epoxy  cannot  be  large.  This  is 
due  to  the  large  ratio  between  axial  and  transverse  Young’s  moduli  for  this 
material.  Results  for  orthogonal  cracking  are  practically  the  same.  The 
situation  is  very  different  for  the  shear  modulus.  Results  for  orthogonal 
cracking  are  unfortunately  not  available  in  this  case. 

We  now  turn  to  numerical  methods  for  cracked  laminate  analysis.  It  is 
evidently  impractical  to  consider  the  case  when  cracks  are  not  equidistant  for 
this  would  require  an  enormous  computer  program.  It  is  for  this  reason  that 
all  numerical  methods  for  composites  are  based  on  idealized  periodic  models. 
For  example,  unidirectional  fiber  composites  have  been  modeled  as  periodic 
hexagonal  arrays  of  Identical  circular  fibers.  The  central  task  is  then  to 
formulate  the  boundary  conditions  on  the  surface  oT  the  repeating  element  and 
this  can  only  be  done  from  symmetry  considerations. 

Consider  for  example  the  orthogonally  cracked  laminate  of  Fig.  7.  In  the 
event  that  the  two  crack  families  are  equidistant  the  repeating  element  is  as 


shown  in  Fig.  8.  When  such  a  laminate  is  loaded  by  constant  biaxial  membrane 
forces  all  the  faces  of  the  repeating  element  and  also  its  midplane3  are 
planes  of  symmetry.  It  is  convenient  for  our  purpose  to  choo3e  as  repeating 
element  the  region  OSx£a;OSySb;OSzSh.  It  follows  that  the 
shear  stresses  must  vanish  on  all  faces,  providing  two  boundary  conditions  on 
each  face.  The  third  condition  for  each  face  is  as  follows:  planes  x,y,z«0, 
normal  displacement  component  vanishes;  plane  z»h,  normal  stress  vanishes; 

planes  x«a,  y-b,  normal  stress  vanishes  on  cracked  part,  normal  components  of 

\ 

displacement  constant,  equal  6a ,  62,  say,  on  uncracked  parts. 

Now  define  the  two  boundary  problems 

u^Ca.y.z)  -  1  t  £  z  S  h  ux(a,y,z)  -  0  t  S  z  S  h 

(1H) 

Uy(x  ,b,z)  -  0  0  S  z  S  t  uy(x,b,z)  -  1  0  £  z  £  h 

while  all  other  listed  boundary  conditions  remain  the  same.  The  elasticity 
solutions  of  these  problems,  satisfying  all  required  traction  and  displacement 
continuity  conditions  on  the  interface  z»t,  may  be  found  numerically  and  are 
labeled  by  superscripts  I  and  II.  It  follows  from  linearity  that  the  fields 
in  the  repeating  element  are  given  by 

u  -  6,  tj1  +  fijU11  o  -  SiO1  ♦  SjO11  (15) 

*% 

In  order  to  determine  the  quantities  5,,62  consider  the  loading  a  ,  o  on  the 

xx  y  y 

laminate  edges.  Then  it  is  easily  shown  that 
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where  overbars  on  the  right  side  denote  averages  over  the  repeating  element 
which  are  known  from  the  numerical  solution.  This  completes  the  field 


solution  and  effective  properties  can  now  be  determined  by  relation  of  stress 

\ 

to  strain  average?  of  repeating  element. 

Note  that  this  procedure  is  by  no  means  restricted  to  cross-plies.  Fig. 

1*1  shows  a  repeating  element  for  a  [0®/0®/-0°]  laminate  which  is  blaxially 

m  n  ns 

loaded  and  which  has  symmetric  interlaminar  crack  distributions  within  the  ±  0® 
plies.  This  can  be  formulated  by  similar  symmetry  arguments  and  the  basic 
solutions  I  and  II  must  now  satisfy  zero  traction  boundary  conditions  on 
internal  crack  surfaces.  But  the  problem  posed  is  very  complex  even  from  a 
numerical  point  of  view  since  the  material  axes  of  the  plies  and  the  cracks 
are  now  at  different  non- orthogonal  directions.  Note  that  the  finite  elements 
must  be  suitable  to  model  the  edges  and  the  cracks  and  in  addition  must  have 
congruent  faces  across  ply  Interfaces  because  of  required  satisfaction  of 
continuity  conditions  at  these  locations. 

Numerical  work  which  has  appeared  in  the  literature  has  been  restricted 
to  singly  and  orthogonally  cracked  cross-plies.  Altus  and  Ishai  [2]  have 
presented  a  finite  difference  analysis  of  singly  cracked  cross-ply  and  have 
given  results  for  in  plane  -  normal  stress  only.  Similar  analysis  and  results 
have  been  obtained  by  Wang  [31]  by  use  of  finite  elements.  Highsmith  and 
Reifsnider  [1H]  have  performed  a  numerical  analysis  of  orthogonally  cracked 
cross-ply  under  tension  on  the  basis  of  a  higher  order  laminate  theory  which 
Is  due  to  Pagano  [20],  which  assumes  linear  variation  of  stress  through  ply 
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thickness.  This  is  not  unlike  the  variational  method  which  is  also  based  on 
gentle  through-thickness  stress  variation,  which  in  view  of  the  two  methods 
are  in  many  respects  similar,  but  the  computational  effort  required  by  the 
numerical  work  is  higher  by  order  of  magnitude. 

3.2.  Thermal  Effects 

Temperature  changes  are  a  major  cause  of  interlaminar  damage.  Consider 
for  example  the  case  of  temperature  change  of  a  [0o/90o]s  laminate.  Cooling 
of  the  laminate  produces  compressive  stresses  in  fiber  direction  and  tensile 
stresses  transverse  to  the  fibers,  of  equal  magnitude,  in  the  plies.  Since 
transverse  tensile  strength  is  much  lower  than  axial  compressive  strength  the 
first  damage  to  appear  is  intralaminar  cracking  of  all  plies,  thus  orthogonal 
cracking.  We  note  in  passing  that  all  laminates  which  are  thermally 
isotropic,  thus  [0°/90°]  [0°/90°/«»5V-*»5,]„,  [0°/60o/-60°]  ,  develop  the  same 

thermal  stresses  and  therefore  the  same  interlaminar  cracking  of  all  plies. 
It  is  important  to  realize  that  when  a  laminate  i3  heated  the  signs  of  ply 
stresses  reverse;  the  fiber  direction  stress  is  tensile  and  the  transverse 
stress  is  compressive.  Therefore  heating  does  not  produce  interlaminar 
cracking. 

The  effective  thermal  expansion  coefficients  (ETEC)  a*j  of  any 
heterogeneous  body  are  defined  by  the  relation 

'u  '  “Ij 8  m 

where  are  the  average  strains  produced  by  the  temperature  change  9.  Note 
that  (17)  merely  states  that  the  average  strains  are  proportional  to  the 
temperature  change.  It  follows  at  once  that  the  ETEC  of  a  cracked  laminate 
may  be  determined  by  evaluation  of  the  average  strains  due  to  temperature 
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change  In  such  a  laminate.  This  is  a  very  difficult  problem  of 
thermo  elasticity  and  it  can  fortunately  be  avoided. 

To  see  this  consider  a  cracked  laminate  which  is  subjected  to  constant 
membrane  load  and  to  temperature  change.  In  this  case  the  average  strains  are 

the  sum  of  (1)  and  (17).  Denoting  the  internal  stresses  due  to  load  alone  as 

s  t 

and  the  stresses  due  to  unit  temperature  change  alone  by  we  have  from 

superposition 
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Levin  [18]  has  proved  a  remarkable  theorem  which  states 


Jv  °u(!)dV  ’  V  (19) 

It  is  seen  that  the  ETEC  are  determined  in  terms  of  the  known  constituent 
thermal  expansion  coefficients  and  the  stresses  due  to  mechanical  load.  The 
remarkable  aspect  of  the  theorem  is  in  that  the  thermal  stresses  do  not  enter. 
Now  the  mechanical  stresses  have  been  discussed  before  and  it  is  thus  seen 
that  stress  analysis  of  a  loaded  cracked  laminate,  whether  analytical  or 
numerical,  also  provides  the  ETEC. 

The  relation  (19)  has  been  exploited  in  [11]  to  obtain  in  simple  fashion 

r% 
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the  ETEC  of  an  orthogonally  cracked  cross-ply,  using  for  the  approximate 

variational  stresses  obtained  in  [10].  Fig.  15  shows  the  variation  of  ETEC 
with  crack  density  for  T300/Epoxy  plies,  assuming  equal  crack  density  in  all 


plies.  Also  shown  are  results  when  only  one  kind  of  plies  Is  cracked.  It  is 
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seen  that  ETEC  diminish  substantially  with  increasing  crack  density  and  that 
they  approach  asymptotically,  the  fiber  direction  TEC  of  the  ply.  This  last 
result  has  been  proved  analytically  [11]  and  its  physical  interpretation  is  as 
follows:  for  large  crack  density  the  plies  lose  their  transverse  stiffness 

(E.J.-0)  and  therefore  the  90°  plies  no  longer  affect  the  axial  expansion  of  the 
0°  plies. 

Loading  and/or  temperature  change  produce  either  tensile  or  compressive 
transverse  stresses  in  the  plies.  The  former  open  cracks  while  the  letter 
close  existing  cracks.  This  produces  complex  thermo- el  as  tic  behavior  since 
there  are  several  possibilities:  all  cracks  open,  some  cracks  open  and  all 
cracks  closed.  In  the  last  case  the  ETEC  are  simply  those  of  the  uncracked 
laminate.  Fig.  16  shows  schematically  possible  values  of  ETEC  for  combination 
of  uniaxial  load  (tension  or  compression)  and  temperature  change  (heating  or 
cooling),  [11]. 

Adams  and  Herakovich  [1]  and  Herakovich  and  Hyer  [12]  have  performed 
numerical  analysis  based  on  finite  elements  to  obtain  ETEC  of  cross-plies  in 
which  the  90°  plies  are  cracked.  They  also  obtained  significant  reduction  of 
ETEC  with  increasing  crack  density. 

3.3.  Damage  Prediction.  Failure  Analysis. 

So  far  we  have  considered  the  problem  of  the  relation  of  properties  to 
damage  in  the  form  of  crack  density.  Since  crack  density  is  not  easily 
measured  it  is  natural  to  inquire  if  it  can  be  analytically  predicted.  Initial 
attempts  to  do  this  [24,25]  were  based  on  the  assumption  that  cracks  will 
form  saturation  patterns  (characteristic  damage  state)  at  intercrack  distances 
which  are  large  enough  to  neglect  crack  interaction.  But  later  analytical  and 
experimental  evidence  does  not  seem  to  support  this  assumption.  The  analysis 
given  in  [9]  for  the  singly  cracked  cross-ply,  while  approximate,  is  certainly 
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much  more  accurate  then  the  simple  shear  lag  approach.  According  to  this 
analysis  crack  interaction  effects  can  be  disregarded  when  the  ratio  of  crack 
interdistance  to  90°  ply  thickness,  p  say,  is  larger  than  M .  But  according  to 
experiments  reported  in  [13]  for  glass/ epoxy  cross-plies  [0°/90°]a  and 
[0°/90j]s  the  value  of  p  is  as  low  as  1.1  and  is  between  the  values  1.1  and  1.6 
for  substantial  ranges  of  the  intralaminar  cracking  process.  But  at  such  low 
values  of  p  interaction  between  cracks  is  very  significant. 

Very  simple  analytical  considerations  show  that  the  highest  in-plane 
stress  oxx  in  load  direction  in  the  cracked  90°  ply  is  midway  between  cracks. 
It  would  thus  be  tempting  to  create  a  very  simple  model  of  crack-density  to 
load  relation  according  to  which,  with  the  help  of  the  analysis  of  [9], 
formation  of  new  cracks  would  occur  when  the  load  is  large  enough  for  the 
mid-crack  stress  to  attain  the  transverse  tensile  failure  stress  of  the  ply. 
But  the  realism  of  such  a  procedure  is  doubtful  since  the  failure  stress  is 
generally  not  some  definite  number  but  a  random  variable  the  reason  Tor  that 
being  that  formation  of  a  crack  is  due  to  a  local  micro-defect  and  is  thus  a 
random  event  whose  probability  increases  with  increasing  stress.  Therefore 
probability  of  occurrence  of  a  new  crack  in  the  region  between  two  existing 
ones  depends  not  only  on  the  value  of  maximum  stress  but  also  on  its  variation 
in  the  interval.  Analysis  requires  assumption  of  some  form  of  probability 
function.  An  example  of  such  work  is  given  in  [2]  which  also  cites  previous 
work  on  thi3  subject. 

Everything  discussed  so  far  is  based  on  the  simplification  of  constant 
through-thickness  in-plane  stress  or  in  other  words:  the  important  quantity  in 
a  thin  ply  is  the  in-plane  force.  If  the  ply  is  thick,  stress  variation  through 
the  thickness  can  become  Important  and  a  more  refined  approach  may  be  needed. 
Wang  [30]  has  considered  this  problem  in  terms  of  critical  flaws  in  the  form 
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of  transverse  micro-cracks  whose  lengths  may  be  much  smaller  than  ply 
thickness.  This  approach  requires  intricate  numerical  analysis  combining  crack 
criticality  and  probability  concepts.  Monte  Carlo  simulation  of  intralaminar 
crack  density  in  singly  cracked  graphite/epoxy  cros3-plies,  produced  by  static 
load  was  found  to  be  in  very  good  to  reasonable  agreement  with  experimentally 
determined  crack  density.  But  we  should  remember  the  thin  ply  problem 
difficulties  discussed  above.  The  Implication  for  the  present  subject  is  that 
a  transverse  crack  length  must  be  much  larger  than  UFC  RVE  in  order  sto  be 
amenable  to  fracture  mechanics  treatment  and  it  would  therefore  seem  that 
this  approach  should  be  confined  to  plies  which  are  sufficiently  thick. 

It  is  at  present  an  open  question  whether  damage  prediction  methods  can 
be  carried  beyond  the  simple  case  of  single  ply  cracking  and  whether  they  can 
develop  into  an  engineering  tool. 

The  problem  of  damage  prediction  becomes  even  more  difficult  with  crack 
saturation  and  the  development  of  interlaminar  damage.  If  one  adopts  the 
point  of  view  that  the  opening  of  interlaminar  cracks  is  a  a  problem  which  is 
governed  by  the  laws  of  classical  fracture  mechanics,  then  it  becomes 
necessary  to  formulate  a  criterion  for  interlaminar  crack  criticality  and  to 
employ  it  to  predict  branching  of  intralaminar  cracks  into  the  interface.  If 
there  are  Intralaminar  cracks  in  neighboring  plies,  e.g.  orthogonal  cracks  in  a 
cross  ply,  the  crossing  points  of  cracks  in  the  different  plies  become  sources 
of  intralaminar  damage.  A  numerical  finite  element  procedure  for  such 
delamination  problems  has  been  reported  in  [31]. 

We  should  remember,  that  as  discussed  above,  finite  elements  can  not  be 
as  small  a3  'we  please;  they  must  be  large  enough  to  represent  the  ply 
properties  and  this  is  a  problematic  limitation  on  such  numerical  analyses. 
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Given  the  severe  complexities  of  damage  and  failure  prediction  in  terms 
of  micro-analysis  and  adding  the  fundamental  difficulty  of  finite  element  3lze 
limitation  it  would  appear  to  us  that  an  engineering  tool  for  damage  and 
failure  prediction  could  hardly  be  formulated  on  the  basis  of  MIDM  alone  and, 
as  we  have  said  in  the  beginning,  a  complementary  combined  MIDM  -  MADM  is 
needed.  A  simple  example  of  such  an  approach  i3  provided  by  the  work  of 

Poursartip  et  al.  [22]  who  have  developed  a  method  to  evaluate  the  fatigue 

\ 

life  of  a  laminate  by  description  of  damage  evolution  as  a  continuum  damage 
process,  thus  MADM,  where  the  necessary  ingredient  to  perform  the  analysis  is 
the  stiffness  redaction  of  the  laminate,  which  can  be  obtained  by  MIDM 
methods.  Further  interesting  developments  in  such  directions,  within  the 
context  of  fatigue  of  laminates,  can  be  found  in  [19,23]. 

<4.  Conclusion 

Most  of  our  discussion  has  been  concerned  with  damage  mechanics  of  fiber 
composite  laminates,  a  subject  of  significant  engineering  importance.  Most  of 
the  work  discussed  has  been  of  MIDM  nature  since  most  of  the  literature  on 
laminate  damage  has  been  based  on  this  approach.  It  may  be  mentioned  that,  in 
contrast,  most  of  the  work  on  metal  damage  has  followed  the  MADM  or 
continuum  approach.  The  reason,  we  believe,  is  in  that  the  internal  damage 
geometry  of  laminates  is  very  much  simpler  than  in  metals,  where  one  ha3  to 
deal  with  irregular  grain  boundaries. 

Even  so  the  micro-analysis  problems  which  must  be  solved  are  very 
complex,  and  assume  Increasing  complexity  with  advancing  stage  of  damage  and 
with  engineering  realism  of  loading  patterns  and  laminate  stacking  sequence. 
It  is  therefore  not  surprising  that  most  of  the  literature  has  been  concerned 
with  analysis  of  symmetric  cross-plies,  the  simplest  laminates  In  existence. 


Apart  from  the  mathematical  difficulties  arising  in  problems  which 
involve  many  Interacting  cracks  there  are  difficulties  of  even  more 
fundamental  nature.  These  have  been  brought  out  in  our  discussion  of  the  thin 
ply  problem  and  the  associated  difficulty  of  smallest  permissible  size  of 
finite  elements  representing  a  heterogeneous  medium.  And  there  is  also  the 
question  to  what  extent  are  micro-cracks  in  a  composite  amenable  to  treatment 
by  classical  fracture  mechanics  methods  in  terms  of  the  effective  elastic 
properties  of  same  composite.  v 

We  would,  in  conclusion,  venture  the  opinion  that  MIDM  for  damaged 
laminates  should  be  pursued  to  the  best  of  our  ability,  that  MADM  should  be 
developed  in  concert  with  MIDM  in  order  to  deal  with  those  cases  which  are  too 
complicated  for  MIDM  treatment  and  that  efforts  should  be  made  to  develop 
methods  to  deal  with  the  fundamental  difficulties  outlined  above. 
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Fig.  9  -  Stiffness  reduction  of  (0°/90°]s  glass/epoxy  laminate.:  theory  and  experiment 


Variation  of  Poisson's  ratio  with  crack  density;  singly  and  orthogonally  cracked  graphite/epoxy  [0°/90°] 
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THERMAL  EXPANSION  COEFFICIENTS  OF  CRACKED  LAMINATES 
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ABSTRACT 

Effective  thermal  expansion  coefficients  of  cross-ply  laminates  containing 
many  intra-laminar  ply  cracks  are  evaluated  as  functions  of  ply  thermo¬ 
mechanical  properties  and  crack  density.  The  method  employed  is  applicable  to 
any  cracked  laminate  if  internal  stress  fields  due  to  mechanical  loads  are 
known.  Results  obtained  demonstrate  strong  reduction  of  expansion 

coefficients  due  to  cracks.  It  is  shown  that  because  of  opening  of  cracks  by 
tension  and  closing  by  compression  the  laminate  has  several  different 
expansion  coefficients,  each  valid  for  a  different  load-temperature 


combination. 
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1.  INTRODUCTION 

There  has  been  considerable  interest  in  recent  years  in  the  analysis  of 
property  deterioration  of  damaged  laminates.  The  primary  damage  consists  of 
many  intra-laminar  ply  cracks,  in  fiber  directions,  and  the  m^jor  research 
effort  has  been  devoted  to  evaluation  of  stiffness  reduction.  Reifsnider  and 
coworkers  have  applied  simple  shear  lag  methods  (see  e.g.  Highsmith  and 

Reifsnider  [1]);  Laws  and  Dvorak  [2]  have  used  so-called  self  consistent 
approximations  to  estimate  cracked  ply  properties  and  on  that  basis  the 

cracked  laminate  stiffnes  Highsmith  and  Reifsnider  [3]  gave  a  numerical 
analysis  of  internal  stresses  and  stiffness  of  cross-ply  laminates  in  which 
both  plies  are  cracked,  based  on  higher  order  laminate  theory.  The  same 

problem  has  been  treated  by  Wang  [U]  with  finite  elements.  Hashin  [5,6]  used 
variational  methods  to  analyze  stiffness  reduction  and  internal  stresses  in 
cross-plies  cracked  in  one  or  in  two  directions. 

By  contrast  the  practically  very  important  problem  of  the  modification  of 
thermal  expansion  coefficients  (TF.C)  by  intra-laminar  cracks  has  received  only 
scant  attention.  Indeed  the  only  published  work  on  this  subject  known  to  us  is 
a  finite  element  numerical  analysis  by  Herakovich  and  Hyer  [7]  who  gave 

results  for  a  cross-ply  cracked  in  two  directions. 

In  the  present  work  we  show  that  in  order  to  evaluate  the  TEC  of  any 
cracked  laminate  it  is  sufficient  to  know  the  internal  stresses  produced  by 
mechanical  load.  Thus  the  TEC  of  a  cros3-ply  cracked  in  two  directions  are 
evaluated  in  simple  fashion  on  the  basis  of  the  stresses  obtained  in  ref.  [6]. 

Since  temperature  changes  produce  both  tensile  and  compressive  internal 
ply  3tresse3  in  laminates  and  since  compression  closes  cracks,  the  effective 
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TEC  of  a  cracked  laminate  will  have  different  values  depending  on  whether  all 
cracks  are  open,  some  are  closed  or  all  are  closed.  This  complex  situation 
will  be  illustrated  for  the  cross-ply  laminates  to  be  analyzed. 


2.  GENERAL  DEVELOPMENT 


To  derive  the  effective  TEC  of  a  composite  material  let  a  body  made  of 
such  material  be  subjected  to  the  homogeneous  traction  boundary  conditions 


Ti(S)  -  o"n, 


where  a.  ■  is  a  constant  stress  tensor  and  n.  are  the  components  of  the  outward 
^  J  J 

normal,  and  to  a  uniform  temperature  change  0  relative  to  a  uniform  reference 
temperatu-e  T  .  Then  from  the  average  stress  theorem  (see  e.g.  Hashin  [8]) 


o 

0.  .  •  <j.  . 

ij  U 


where  over  bar  denotes  volume  average.  Then  the  average  strain  is  given  by 


hi  ■  s‘jkiV  ‘  "ij9 


where  S?.,,  is  the  effective  compliance  tensor  and  af.  is  the  effective  thermal 
ijkl  ij 

expansion  tensor. 

Let  the  stresses  due  to  load  (2.1)  only,  with  no  temperature  change,  be 
denoted  and  let  the  stresses  due  to  unit  uniform  temperature  rise  only, 

with  zero  boundary  tractions  T^(S)*0,  be  denoted  oh-  If  follows  by 
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s  iper  posi  tion  that  the  stress  field  in  the  composite  material  due  to  load  and 
temper  at  ire  change  can  be  expressed  in  the  form 


Oij(x)  -  ofjW  +  0  0^00 


Levin  [9]  has  proved  a  remarkable  theorem  which  states 


For  a  multiphase  body  with  homogeneous  constituents  (2.5)  assunes  the  form 


Y  (m)-(m)  *  o 

^“ij  °  ij  Vm  "  aij°ij 


where  m  is  the  phase  index,  v^  the  volume  fraction  of  mth  phase  and 
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The  effective  TEC  can  be  explicitly  expressed  from  (2.6)  in  the  form 


*  y  (m)  (m) 
an  =  c_a  i  B .  ,  .  .  v 
ij  k1  klij  m 


where  B  ™  is  defined  by  the  linearity  relation 
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The  remarkable  aspect  of  the  relations  (2.5),  (2.6),  (2.9)  Is  In  that  they 

g 

depend  only  on  the  mechanical  stresses  o.^  and  not  on  the  thermal  stresses 


The  relation  (2.9)  has  been  exploited  by  Levin  (1967)  to  derive  direct 
relations  between  the  effective  elastic  properties  and  effective  TEC  of  twov 
phas e  materials.  In  the  case  of  composite  materials  which  consist  of  more 
than  two  phases  such  direct  relations  cannot  be  established,  but  evaluation  of 
effective  elastic  froperties  is  still  sufficient  for  evaluation  of  effective 
TEC  for  the  former  requires  determination  of  o®.(x)  which  from  (2.5)  and  its 
corollaries  is  sufficient  to  determine  <j*  . 

We  shall  now  consider  the  application  of  these  relations  to  evaluation  of 
TEC  of  damaged  laminates.  For  sake  of  simplicity  we  consider  only  symmetric 
laminates  which  are  thermoelas tically  orthotropic.  Generalization  to  more 
complicated  cases  is  a  str aightf orw ard  matter. 

A  laminate  of  rectangular  plan  form  with  orthotropy  axes  parallel  to  the 
sides  of  the  rectangle  is  sutyected  to  membrane  load3 


Nxx  "  2h0 


N  -  N 

yy  xy 


(2. 10) 


where  2h  is  laminate  thickness,  fig.  1.  The  present  loading  is  a  special  case 
of  (2.1)  when  the  only  surviving  component  of  o°j  is  o°* =  0^  1 . 


Let  the  internal  stress  field  in  the  damaged  laminate  be  represented  as 


Let  the  internal  stress  field  in  the  damaged  laminate  be  represented  as 


ofjOO  -  0  ij(— )  v 


(2.  11) 


where  o^tx)  is  the  stress  field  in  the  undamaged  laminate  and  ojj(x)  is  the 
change  in  stress  due  to  the  internal  defects.  Introducing  (2.11)  into  (2.5)  we 
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(2.12) 


Since  (2.5)  is  also  valid  for  the  undamaged  laminate  we  identify  the  first 
integral  in  (2.12)  as  the  TEC  a”x  of  the  undamaged  laminate. 

It  is  well  known  that  the  a  ^  are  constant  in  each  layer  or  ply.  Let 
these  stresses  in  the  mth  layer  be  referred  to  the  material  axes  of  the 
unidirectional ly  reinforced  ply  material  and  be  expressed  in  the  form 


(m)  o  (m) 

0  y  ’  0  k  u 


(2. 1  3) 


It  follows  at  once  that 


E(m).  (m) 

0  u k  u 


(2.  I  4) 


and  it  is  readily  realized  that  k  are  components  of  the  B  tensor  in  (2.9) 
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for  the  present  case. 

In  a  more  general  sense  we  conclude  that  for  a  damaged  laminate,  and 
indeed  for  any  damaged  composite  material 


*  o 

a.  .0 .  , 
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o  o 
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(2. 15) 


where  is  the  stress  perturbation  due  to  defects.  v 

We  arrive  at  an  Immediate  interesting  consequence  by  noting  the  well 
known  fact  that  interior  interlaminar  cracks  have  no  effect  on  the  stress 
distribution  in  symmetric  laminates  under  membrane  loading  since  such  cracks 
have  no  effect  on  the  states  of  plane  stress  in  the  layers.  It  follows  that 
a^.-0  and  therefore  from  (2.15)  30(1  thus  the  effective  TEC  in  the 

laminate  plane  are  not  affected  by  interlaminar  cracks. 

3.  EVALUATION  OF  THERMAL  EXPANSION  COEFFICIENTS 
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We  shall  here  be  concerned  with  the  case  of  a  symmetric  cross-ply 
laminate,  denoted  [0/90°]  in  usual  laminate  notation  in  which  one  kind  (e.g. 
90°)  or  both  kinds  of  plies  contain  distributions  on  intra-lamlnar  cracks,  fig. 
1.  The  case  of  cracked  90°  occurs  when  the  laminate  is  loaded  normal  to  the 
90  layer.  The  case  of  both  plies  cracked  occurs  when  the  laminate  is 
subjected  to  he  a  ting/  cooling  temperature  cycles,  for  in  such  a  case  the 
stresses  normal  to  the  fibers  in  the  different  plies  are  equal  and  opposite  in 
sign.  Since  cracks  are  produced  by  tensile  stresses  heating  can  produce  cracks 
in  one  ply  and  cooling  can  produce  crack3  in  the  other. 


l 
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Another  way  to  produce  orthogonal  cracks  in  both  plies  is  to  subject  a 
±45°  laminate  to  cyclic  tension-tension  load  in  direction  of  the  bisector  of 
the  ±45°  directions. 

g 

To  evaluate  effective  TEC  we  need  the  o^j  stress  field  in  an  orthogonally 
cracked  cross-ply  laminate.  This  problem  has  been  treated  on  the  basis  of 
finite  elements  by  Wang  [4]  and  on  the  basis  of  a  higher  order  laminate  theory 
by  Highsmith  and  Reif snider  [3],  both  of  these  treatments  being  numerical.  At\ 
analytical  approximate  variational  treatment  of  the  problem  has  been  given  by 
Hashin  [6]  and  we  shall  use  these  results  to  evaluate  effective  expansion 
coefficients  analytically.  It  is  of  interest  to  note  that  the  stress  field 
used  is  admissible  in  the  sense  of  the  principle  of  minimum  complementary 
energy  and  it  can  thus  be  used  to  obtain  a  strict  lower  bound  on  effective 
Young's  modulus  of  the  cracked  laminate,  ref.  [6].  The  present  use  of  this 
stress  field  does  not  necessarily  produce  a  bound  on  the  effective  TEC.  To 
establish  such  bounds  the  thermal  stress  field  Oy  is  also  needed  and  this 
will  be  discussed  elsewhere.  The  justification  for  the  present  analysis  is 
outstanding  agreement  of  previous  results  for  Young's  modulus  of  a  cracked 
laminate  with  experimental  results,  [5]. 

It  is  assumed  for  the  sake  of  simplicity  that  the  cracks  in  both  plies  are 
equidistant.  In  the  following  the  inner  ply  is  assigned  the  90°  direction  and 
is  labelled  by  index  1  while  the  outer  plies  are  assigned  the  0°  direction  and 
are  labelled  by  index  2.  The  intercrack  distances  in  the  1  and  2  plies  are 
denoted  2a  and  2b  respectively.  If  the  load  (2.10)  Is  applied  to  the  laminate 
it  follows  from  symmetry  that  the  stress  field  is  periodic.  The  repeating 
element  for  the  stress  field  may  be  chosen  as  the  region  of  the  laminate 
defined  by  4  plane  cuts  along  any  two  intersecting  crack  pairs  and  the  two 


planes  z=0,  h,  fig.  2. 


The  perturbation  stresses  due  to  cracks  as  defined  in  (2.11)  are  denoted 


ofP  and  o  j in  the  regions  1  and  2  respectively.  The  ply  materials  are 


transversely  isotropic  with  TEC  a  in  fiber  direction  and  a_  transverse  to  the 

A  1 


fibers.  It  follows  that  the  last  integral  of  (2.15)  can  be  expressed  in  the 
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According  to  the  analysis  given  in  ref.  [6]  the  stresses  entering  into  (3.1) 


have  the  form 
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^  and  are  the  stresses  in  the  1  ply  of  the  undamaged  laminate 


where  o  and  a 
x  y 


and  the  functions  <j>  and  i}>  are  defined  by  linear  integro- differential  equations, 
[6],  and  the  boundary  conditions 


Since 
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These  integrals  vanish  in  view  of  (3. 3h)  and  therefore  the  ozz  stresses  do  not 
contribute  to  (3.1).  Defining  the  mean  values 
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It  then  follows  in  straightf orward  fashion  from  (3.1)  and  (2.12)  that 
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The  values  of  <j>  and  ^  have  been  determined  in  ref.  [6]  and  are  given  in  the 
appendix. 
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We  now  consider  some  special  cases.  The  first  is  a  laminate  in  which 
only  one  kind  of  plies  is  cracked.  Let  it  be  assumed  that  there  are  no  cracks 
in  the  outer  plies.  This  implies  that  intercrack  distance  b  is  infinitely 
large,  thus  by  definition  ( A— 2)  in  the  appendi  x/  It  follows  from  (A-7)  and 

(A-9)  that  in  this  case  w^O,  then  from  (A-10) 


l  -m^m  1 1  -<*>1  j 


-  n>20(1) 


and  these  values  have  to  be  introduced  into  (3-5)  to  obtain  the  TEC. 

Similarly,  if  the  inner  ply  is  not  cracked  and  the  outer  ones  are,  then 


1  -m  ( 1  -uj^) 


Then  from  (3,5)  and  (A-5)  the  TEC  in  the  two  cases  are 
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(a) 

(3.8) 

(b) 

where  the  constants  AQ>  Bq  and  CQ  are  defined  in  (A-3). 

The  TEC  assume  asymptotic  limit  values  for  large  crack  density.^ 
Theoretically  this  implies  that  either  one  or  both  of  the  ratios  p  ,  in 
(A-2)  goes  to  zero.  This  does  not  really  happen  in  practice  since  the  crack 
patterns  generally  reach  saturation  for  some  definite  intra-crack  distances. 
Still  the  limiting  values  are  of  interest  as  lower  bounds  and  it  will  also  be 
seen  that  they  are  quite  quickly  approached  as  the  crack  density  increases. 

It  is  easily  shown  that  when  p+0,  the  quantity  u>  as  defined  by  (A-7)  and 
by  (A- 9)  approaches  the  limit  1.  It  then  follows  from  ( A- 1 0)  that 


-  12 


*0)  o  “a  aT  .  /r  ..  (1)-(1) 

l  «  O  +  -z ; —  ( 1  +  B  /  C  )k  d> 

XX  XX  1  +1  o  o  x  v 


•a  *C2)  -a°  .  °J^L  (,  .  B  /,  *'■>' 

xx  xx  1  o  o  y  * 


(1 )-(2) 

y 


lim  <j>  -  1 

>■ 

lim  ij)  »  1 


(3.9) 


k 

t 

7Z- 


5* 

ft 


Rather  than  inserting  (3-9)  into  (3.5)  is  is  more  convenient  to  use  (2.12)  with 
(3.9).  It  then  easily  follows  that 
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In  the  event  when  only  the  inner  or  outer  plies  are  cracked  we  have  from 


(3.  8) 
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(3. 1  D 


Evaluation  of  the  TEC  a  and  a ^  is  of  course  quite  similar. 

For  this  purpose  the  TEC  with  subscripts  xx  in  (3.5)  and  (3.8)  should  be 

replaced  with  the  corresponding  TEC  with  subscripts  yy  and  k  ^  and  k  ^ 

x  y 

should  be  interpreted  as  (3- 5b)  when  the  loading  (2.10)  is  replaced  by 

%y-2h°°  <3-9) 

Figs.  3,i4  show  the  variation  of  TEC  a*x  with  crack  density  for  Glass/Epoxy 
and  for  Graphite/Epoxy)  (T300  fibers)  cross-plies.  In  each  case  t^-t^,  X=1  and 
results  are  given  for  orthogonal  cracking,  a=b,  and  also  when  only  the  inner 
layer  cracked.  Ply  elastic  properties  and  stress  coef f icients  are  given  in 
table  1.  It  is  seen  that  the  cracks  have  significant  effect  on  the  TEC  of  the 
laminate,  especially  in  the  case  of  T300/Epoxy;  that  for  orthogonal  cracking 
the  TEC  are  lower  than  for  only  one  ply  cracked  and  that  for  a=b»t1  TEC 
values  are  already  very  close  to  large  crack  density  asymptotic  limits. 

When  a  cracked  laminate  is  subjected  to  loading  and  temper  at  ire  change 
i 

|  the  internal  ply  stresses  can  be  tensile  or  compressive.  Since  tensile 


I 

i 


transverse  (normal  to.  fibers)  stresses  open  the  cracks  and  compressive 
transverse  stresses  close  them,  the  thermal  deformation  of  the  laminate  will 


be  governed  by  abruptly  changing  TEC  during  the  load-temperature  process. 
This  will  now  be  Illustrated  for  a  cross-ply  laminate  in  which  all  plies  are 
cracked. 

The  transverse-stresses  are  and  and  the  signs  of  these 

xx  y  y 

stresses  due  to  load  c°  in  x  direction  or  due  to  temperatire  change,  for 
uncracked  laminates,  can  be  seen  in  table  1.  It  is  probably  safe  to  assume 
that  these  stresses  will  retain  their  signs  (but  of  course  not  their  values) 
when  the  laminate  is  cracked.  Consider  first  the  case  when  a  stress-free  and 
load  free  cracked  cross-ply  is  heated.  In  this  case  the  transverse  stresses  in 
both  plies  are  compressive.  Therefore  the  cracks  close  and  the  TEC  is  that  of 
the  uncracked  laminate.  On  the  other  hand,  when  the  laminate  is  cooled  from 
stress-free  state,  the  transverse  stresses  are  tensile,  the  cracks  open  up  and 
the  TEC  are  given  by  (3*5). 

Next  consider  the  same  laminate  which  is  heated  under  tension  Nxx*2o°h. 
The  load  alone  produces  tensile  transverse  stresses.  Therefore  in  the  initial 
heating  stage  the  TEC  in  x  direction  is  a*x-  At  some  temperature  the 
compressive  thermal  stresses  will  close  the  cracks  in  one  ply,  ply  2  say.  At 

*(  i ) 

this  stage  the  TEC  becomes  axx  ,  eq.  (3,8a).  Further  heating  will  increase 
the  thermal  compressive  transverse  stresses  to  the  point  where  the  cracks  in 
ply  1  will  also  close  and  from  thereon  the  TEC  is  axx,  For  cooling  the  TEC  is 
always  a*x  since  the  thermal  transverse  stresses  are  tensile  and  are  added  to 
the  tensile  mechanical  stresses. 

Similar  reasoning  can  be  applied  to  the  case  of  temperature  change  with 
compression.  The  nature  of  the  thermal  deformation  for  this  and  the  previous 
cases  is  schematically  shown  in  fig.  5. 
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4.  DISCUSSION  AND  CONCLUSION 

It  has  been  shown  that  the  effective  TEC  of  cracked  laminates  can  be 
determined  in  simple  fashion  if  the  internal  stress  fields  produced  by  constant 
membrane  loads  are  known.  Such  stress  fields  can  be  utilized  in  this  fashion 
whether  known  analytically  or  numerically.  In  this  work  we  have  utilized 
simple  approximate  stress  fields,  obtained  elsewhere. 

\ 

<*« 

It  has  been  found  that  intralaminar  cracks  significantly  reduce  the  TEC  of 
the  laminate  and  this  is  in  agreement  with  the  numerical  results  given  in  ref. 
7.  When  both  plies  of  a  cross-ply  laminate  are  cracked  the  limit  of  reduction 
of  TEC  is  the  axial  TEC  of  a  ply.  Therefore  the  TEC  reduction  is  particularly 
significant  for  Graphi te/Car bon  fiber  materials  since  in  these  cases  the  axial 
TEC  of  the  ply  is  very  small  and  often  negative.  This  is  illustrated  in  fig.  4. 

An  interesting  phenomenon  occurs  due  to  the  fact  that  thermal  stresses  in 
the  plies  are  both  tensile  and  compressive.  In  a  cracked  ply  under 
compression  the  cracks  close  and  therefore  such  a  ply  behaves  as  an  un cracked 
ply.  Therefore  the  laminate  TEC  during  a  load-temperatir e  program  will  have 
different  values  and  it  is  necessary  to  evaluate  the  thermal  stresses  to 
determine  which  TEC  is  appropriate,  the  choice  being  TEC  for  both  plies 
cracked,  one  ply  cracked  or  none.  A  variational  analysis  of  thermal  stresses 
will  be  presented  elsewhere. 

Since  a  cross-ply  is  hardly  a  laminate  of  practical  importance  it  is 


EUi 


it 


SB 


appropriate  to  Inquire  into  the  problem  of  TEC  modification  by  cracks  in 
laminates  of  mere  complex  configuration.  On  the  basis  of  present  treatment 
this  problem  is  resolved  when  the  stresses  due  to  membrane  loads  are  known 
and  it  is  therefore  in  this  direction  where  we  have  to  look  for  the  answers. 
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Obviously  balanced  symmetric  0/90  and  +45/_I<5  laminates  with  same  crack 
density  in  all  plies  are  isotropic  in  thermal  expansion  with  same  TEC.  It  is 
therefore  to  be  expected  that  balanced  symmetric  0/90/+H|f/-il5  laminates  with 
same  crack  density  in  all  plies  will  have  these  same  TEC,  at  least 
approximate!  y. 
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Material  properties  of  reinforced  ply 
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-  Axial  Young's  modulus  (fiber  direction) 

-  Axial  Poisson's  ratio 

-  Transverse  Poisson's  ratio 
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Laminate  Stress  Coefficients 
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l(1) 

X 
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.75 
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-.041 1 

2.091 
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PHKflf 

.0733 

.41  31 

-.21  96 

Gi/Ep 

BrH 

.50 

.50 

.4644 

-.0764 
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.0764 
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.1960 

.1960 

-.1960 

.75 

.25 

.21  82 

-.201 

3.345 

.0603 

-.2477 

.0927 

.7430 

-.2783 

T300/Ep 

.50 

.50 

.1221 

-.0325 

1.878 

.0325 

-.2749 

.274  9 

.274  9 

-.2749 

MPa/C° 

Elastic  Moduli 
TEC 


GPa 

10_61  /C° 


Temperature  Stress  Coefficients  MPa/C°  (stress  per  unit  temperature  rise) 
Table  1  -  Thermo-elastic  properties  and  stress  coefficients 
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T300/ EPOXY 


Variation  of  thermal  expansion  coefficient  of  Graphite/Epoxy  cross-ply 

with  crack  density. 
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ABSTRACT 


A  Differential  Scheme  (DS)  approximation  for  elastic  properties  of  cracked  materials  is  esta¬ 
blished  by  a  limiting  process  on  the  basis  of  the  DS  for  porous  materials.  The  method  is  applied  to 
obtain  stiffness  reduction  due  to  randomly  oriented  elliptical  and  penny  shaped  crack  distributions  in 
isotropic  matrix,  and  to  the  case  of  aligned  plane  cracks  in  orthotropic  sheets. 
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I 

INTRODUCTION  I 


The  problem  of  stiffness  reduction  of  elastic  bodies  due  to  the  development  or  presence  of 
many  disordered  cracks  is  of  significant  scientific  and  engineering  importance  and  has  been  the  sub¬ 
ject  of  many  investigations.  This  appears  to  be  a  particularly  difficult  problem  of  the  elastic  beha¬ 
vior  of  heterogeneous  media  as  is  evidenced  by  the  fact  that  the  only  exact  results  derived  to  date 
are  for  the  case  of  small  crack  density,  when  the  cracks  do  not  interact  (Periodic  crack  arrays  can 
be  analyzed  numerically  to  any  desired  degree  of  accuracy.) 

Variational  bounding  methods  which  have  proved  so  useful  for  composite  media  have  been  only 
sparsely  applied  to  this  problem:  Gottesman  (1980)  for  plane  cracks  and  Willis  (1981). 

The  approximate  method  known  as  the  Self  Consistent  Scheme  (SCS)  has  first  been  applied  to 
this  problem  by  Budiansky  and  O'Connell  (1976)  with  special  attention  to  randomly  oriented  ellipti¬ 
cal  and  penny  shaped  cracks  in  isotropic  matrix.  This  has  been  further  developed  by  Horii  and 
Nemat-Nasser  (1983)  to  take  into  account  friction  between  opposite  crack  surfaces.  Hoenig  (1979) 
has  applied  the  SCS  for  aligned  distributions  of  elliptical  cracks  and  Gottesman  et  al.  (1980)  for 
aligned  plane  cracks  in  orthotropic  sheet.  Aboudi  and  Benveniste  (1987)  have  used  a  generalized 
SCS,  whereby  a  cracked  circular  disk  is  imagined  to  be  embedded  in  the  effective  medium,  rather 
than  the  crack  directly,  for  randomly  oriented  plane  cracks.  Benveniste  (1986)  has  also  applied  the 
Mori-Tanaka  approximation  to  cracked  materials.  Other  literature  is  cited  and  discussed  in  the 
above  referenced  papers. 

Another  approximate  method  of  general  nature  is  the  so-called  Differential  Scheme  (DS).  This 
method  is  in  a  certain  sense  related  to  the  SCS  as  will  be  explained  further  below.  To  the  best  of 
our  knowledge  both  methods  originate  with  Bruggeman  (I93S)  in  the  context  of  conductivity  and 
dielectrics  of  heterogeneous  media.  He  named  the  SCS  symmetric  effective  medium  theory,  since  its 
results  are  insensitive  to  phase  interchange,  and  the  DS  -  effective  medium  theory.  The  former  is  a 
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more  fortunate  choice  of  words  that  self  consistent  scheme,  for  self  consistency  merely  implies  that 
the  effective  elastic  moduli  and  effective  elastic  compliance  tensors  as  predicted,  separately,  are 
mutually  reciprocal.  But  this  is  only  a  minimal  necessary  requirement  which  must  be  imposed  on 
any  approximate  theory  and  is  therefore  not  a  suitable  label  of  identification. 

The  purpose  of  the  present  work  is  to  establish  a  general  DS  formulation  for  cracked  materials. 
To  the  best  of  our  knowledge  the  only  DS  treatments  of  cracked  materials  in  the  literature  are  those 
of  Salganik  (1973),  who  applied  the  Roscoe  (19S2)  interpretation  of  the  DS  to  the  case  of  randomly 
oriented  penny  shaped  cracks  in  ad-hoc  fashion,  and  of  Henyey  and  Pomphrey  (1982)  who  essen¬ 
tially  reproduced  the  Salganik  results. 

The  motivation  for  such  a  study  are  some  inherent  problems  with  the  SCS  results.  For  porous 
media  with  spherical  voids  the  SCS  predicts  that  the  effective  elastic  moduli  diminish  linearly  with 
void  volume  fraction  until  they  vanish  for  50%  voids,  Budiansky  (1965),  Wu  (1966).  This  is  theo¬ 
retically  unacceptable  and  is  in  disagreement  with  experimental  results.  It  is  a  direct  consequence  of 
the  phase  symmetry  of  the  SCS  predictions  as  opposed  to  the  considerable  phase  bias  when  one 
phase  is  matrix  and  the  other  inclusions.  Similar  problems  are  encountered  in  the  SCS  results  for 
penny  shaped  cracks,  Budiansky  and  O'Connell  (1976).  The  effective  Young's  modulus  diminishes 
linearly  with  crack  density  parameter  (CDP)  and  the  effective  shear  modulus  nearly  so,  both  moduli 
abruptly  vanishing  for  CDP  of  value  9/16.  This  is  not  acceptable  since  the  moduli  should  vanish 
asymptotically  with  increasing  crack  density.  Furthermore,  linear  variation  with  CDP  is  restricted  to 
small  crack  density  when  the  cracks  do  not  interact  and  can  therefore  not  be  valid  for  large  density 
when  the  cracks  do  interact. 

It  has  been  shown  by  Boucher  (1976)  that  the  DS  approximation  for  porous  media  does  not 
exhibit  the  SCS  problems  and  agrees  quite  well  with  experimental  results.  It  would  therefore  seem 
worthwhile  to  develop  the  DS  method  for  cracked  materials. 
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GENERAL  DEVELOPMENT 

The  SCS  and  the  DS  may  both  be  related  to  some  fundamental  results  in  the  theory  of  elastic 
two  phase  media.  It  is  well  known,  e.g.  Hill  (1963),  that  when  such  a  material  is  subjected  to  aver¬ 
age  strain  <  or  to  average  stress  a  then  the  effective  elastic  moduli  tensor  C*  and  effective  compli¬ 
ance  tensor  S*  can  be  expressed  in  the  form 


C*  -  +  (C*2)  -  C*1*)  Ay¬ 


s'  -  S(1)  +  (S(2)  -  S(1))  Bv 


(a) 


(b) 


(1) 


where  cO).c(2),s(1)  and  S(2)  are  the  moduli  and  compliance  tensors,  respectively,  of  the  phases 
which  have  volume  fractions,  v( ,  v^  and  the  symmetric  fourth  rank  tensors  A  and  B  are  defined 
by  the  linearity  relations 


<-(2)  -  Ae 


o(2)  -  Bo 


(2) 


where  <  ^  and  o^  are  average  strain  and  average  stress  over  phase  2,  the  first  when  t  is  pre¬ 
scribed  and  the  second  when  o  is  prescribed.  Eqs.  (2)  are  equally  valid  when  the  phase  indices  are 
interchanged. 

We  introduce  from  now  on  the  restriction  that  phase  2  is  in  the  form  of  inclusions  which  are 
embedded  in  matrix  1.  Then  (2)  becomes  a  convenient  device  for  resolving  the  case  of  small 
concentration  v^*  c«l,  in  which  case  it  may  be  assumed  that  the  inclusions  are  isolated.  In  this 
case  the  state  of  strain  or  stress  in  any  of  them  is  as  if  it  were  embedded  all  by  itself  in  infinite 
matrix  with  far  field  <  or  o.  Eshelby  has  provided  a  well  known  solution  for  the  general  eilipsoi- 


dal  inclusion  which  amounts  to  determination  of  A  or  B  in  (2)  now  identified  as  A  and  B  for  this 

s  s 

special  case.  Then  the  small  concentration  result  can  be  expressed  in  the  dual  forms. 


C*5  -  cf1*  +  (C*2)  -  A*c 
Sw  -  S(1)  +  (S(2)  -  S(1))  Bsc 


As  «  As(C<2),C(1),g)(a) 
B*  ■  B*(S(2),  S(1),  g)  (b) 


(3) 


where  g  symbolizes  inclusion  geometry. 

In  the  most  commonly  used  version  of  the  SCS  it  is  assumed  that  when  inclusion  volume 
concentration  is  not  small  the  average  state  of  strain  or  stress  in  any  inclusion  may  be  estimated  as  if 
it  were  embedded  by  itself  in  a  homogeneous  medium  to  which  are  assigned  the  effective  elastic 
properties  C*  and  S*  of  the  composite.  The  consequence  of  this  assumption  is 

C*  -C*l)  +(C(2)  -C(,))A*v2 

(4) 

A*  -  AS(cf2),C*,t) 

with  an  analogous  result  for  S* .  There  are  as  many  eqns.  (4)  as  there  are  components  of  C*. 

We  mention  in  passing  that  there  are  other  SCS  versions  in  which  a  composite  element  is 
embedded  in  the  effective  medium.  For  discussion  see  e.g.  Hashin  (1983). 

The  DS  is  based  on  the  notion  of  incremental  construction  of  the  composite  material  by  gradual 
addition  of  infinitesimal  amounts  of  inclusions.  The  fundamental  assumption  is  that  the  effect  of  an 
increment  dv^  on  the  effective  properties  at  current  v2  is  governed  by  the  small  concentration  res¬ 
ults  (3)  in  which  the  matrix  is  assigned  the  effective  properties.  But  the  equivalent  of  the  small 
concentration  c  is  not  simply  dv2  for  since  the  composite  must  be  statistically  homogeneous  at  all 
stages  of  the  construction  process  accomodation  of  new  added  inclusions  requires  rearrangement  of 
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the  current  ones  and  because  of  this  consideration  it  may  be  shown  that  the  quantity  dv^Al-v^ )  is 
the  proper  equivalent  of  small  concentration.  This  procedure  results  in  differential  initial  value 
problems  for  effective  moduli  and  compliances  which  may  be  expressed  in  the  general  tensor  form, 
McLaughlin  (1977), 


3^-  -  (C*2)  -C*)  A*/(1-v2) 

c*  «  <fl> 

v2-0 

-  (S{2)  -S*)B*/(1-v2) 

s*|  -  s<‘> 


(5) 


(6) 


Self-consistency  of  C*  and  S*  as  predicted  by  (5,6)  has  been  stated  by  McLaughlin  (1977).  A  proof 
is  here  given  in  appendix  I. 

The  incremental  procedure  described  above  has  been  initiated  by  Bruggeman  (1935)  and  has 
been  further  developed  by  Roscoe  (1952),  Boucher  (1976)  and  McLaughlin  (1977).  In  Roscoe's  in¬ 
terpretation  of  the  DS  it  is  assumed  that  there  are  many  inclusion  sizes  of  many  orders  of  magni¬ 
tude.  The  filling  process  starts  with  the  smallest  inclusions  and  proceeds  with  increasing  orders  of 
magnitude.  The  mathematical  expression  of  this  procedure  results  again  in  eqns.  (5,6).  The  reason 
for  this  interpretation  is  that  if  the  added  inclusions  are  much  larger  than  the  current  ones  then  the 
fundamental  assumption  that  they  ’see”  an  effective  medium  is  quite  accurate.  But  it  should  be 
noted  that  the  requirements  of  infinitesimality  of  inclusion  volume  increment,  large  number  of  in- 


elusions  and  increasing  order  of  magnitude  of  inclusion  size  are  mutually  contradictory. 

Norris  (1985)  has  pointed  out  that  there  are  many  ways  of  incremental  construction  of  a  parti¬ 
culate  composite,  each  of  which  leads  to  different  effective  elastic  moduli.  The  physical  fact  is  that 
the  effective  elastic  properties  of  a  composite  are  uniquely  determined  by  its  phase  geometry  (in  all 
details)  and  by  the  phase  properties.  If  the  DS  were  an  exact  procedure  then  such  "path  depen¬ 
dence*  would  be  unacceptable.  Since,  however,  the  DS  is  approximate  its  results  are  of  uncertain 
accuracy,  but  it  is  an  open  question  whether  such  uncertainty  can  be  quantified  by  artificial  variety 
in  filling  process.  In  any  event  the  option  of  filling  path  dependence  does  not  apply  in  the  case  of  a 
cracked  solid. 

In  preparation  for  cracks  we  first  consider  the  case  when  the  inclusions  become  voids.  "This 
implies  that  C<2>  -*  0;  S^-»  oo.  Then  (la)  assumes  the  form 


c'  -c"’*!- Vj) 


where  I  is  the  fourth  rank  symmetric  unit  tensor  and  Aq  is  A  of  (2a)  for  voids.  The  interpretation 


-f2) 

of  e  in  (2a)  in  this  case  follows  from  the  average  strain  theorem 


[V  ■  =77-  (u.  n.  +  u.  n.  )ds 

ij  2V  L  i  J  J 


where  is  the  interface,  thus  the  surface  bounding  the  voids. 

Modification  of  (lb)  for  voids  needs  some  manipulation  because  S^-*  oo.  Introducing  (2b) 
into  (lb)  and  remembering  that  -*  0  for  voids,  we  obtain 


S*-S(,)+e'(2)v, 


Ok  of i ofo f  of alkJ 
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which  can  be  written  in  the  form 


S*  -S(1)(I  +  Bqv2) 


where  the  tensor  Bq  is  defined  by 


B  a 
o 


From  (7,10)  there  immediately  follow  the  small  concentration  results 


i 

i 

i 


(9) 


(10) 


C*  -C(1)(I-  Aqc2) 
S*  -S(1)(I  +  Be) 


(ID 


where  A  and  B  are  determined  by  evaluation  of  (8)  for  isolated  voids  for  far  fields  e  or  a,  res- 
o  o 

pectively. 

Cracks  are  defined  as  very  fiat  voids  of  vanishing  thickness  and  thus  also  of  vanishing  volume. 
We  multiply  both  sides  of  (8)  by  v^  and  consider  the  limit  of  flattening  out  into  cracks,  which  will 
from  now  on  be  identified  as  v^  —  0.  Then  the  surface  integral  for  each  crack  is  confined  to  the 
two  congruent  crack  surfaces  and  thus  each  component  of  normal  appears,  twice  with  opposite  signs. 
Therefore 


lint  («^v  ) 

»  2  v^— *  0 


2V 


[  (tuj  lnj  +  [Uj  ]n.  )ds 
JSc 


»J 


(12) 


where  S  denotes  the  crack  surface  (once),  fu.  1  is  the  displacement  jump  across  the  crack  surface, 
c  i 


i 

V 

s 

l 

» 

I 


usually  called  the  crack  opening  displacement  (COD),  and  n.  are  the  components  of  the  normal  with  ; 

respect  to  one  side  of  the  crack. 

The  tensor  7  in  (12)  is  obviously  linearly  related  to  applied  <  or  5,  whichever  the  case.  We 

I 

write  these  linearity  relations  in  the  form 


(a) 

(13) 

(b) 

where  the  alternate  forms  (13b)  will  be  used  according  to  convenience. 

It  follows  from  (13),  the  left  of  (12)  and  from  the  definition  of  Aq  and  Bq  that 


1  -  ?e 

1°  -  S(1)Qo  «  R  a 


P  »  lim  (A  v, ) 

0  1  )  v  -►  0 
Q-lim(Bov2)  1  v2 


(14) 


It  follows  at  once  that  for  the  case  of  cracks  (7,  9)  transform  into 


C*  -  C(,)  (I  -  P) 

S*  -  S(1)(I  +  Q)  -  S(l)  +  R 


(15) 


These  are  general  results  for  cracked  elastic  bodies,  expressing  effective  elastic  properties  in  terms  of 
COD.  The  case  of  small  crack  density  is  defined  by  non  interacting  cracks.  This  is  a  special  case 
of  (IS)  expressed  as 


c”.C<l)(I-P!) 


(a) 


ranwvwii/wrwnii  w  v 
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SW  -  S(1\l  +  Qs)-S(1)  ♦  R* 


P*  -P*(S(1),gJ 


QS  -  QS(S(1),  g  )  Rs  -  S(1)Qs(S(1),g  ) 

v  V 


where  g  is  crack  geometry.  The  tensors  and  Qs  are  determined  by  evaluation  of  (13)  for  non- 
c 

interacting  cracks. 

It  is  easily  realized  that  the  SCS  for  cracked  bodies  can  be  written  in  terms  of  (13,16)  in  the 
alternative  forms 


C*  -C^fl  -  P*) 


S*  -S(1)(I  +  Q*) 


P*  *  P* (C* ,gc )  Q*  -QS(S*.gc) 


We  now  consider  the  DS  and  start  with  the  case  of  voids.  We  can  easily  apply  the  reasoning 
leading  from  (3)  to  (5,6)  for  the  case  of  voids  by  using  small  concentration  versions  of  (11).  Thus 


c”  -  C0)  [l  -  aJc2  ] 


s“  -  S(,)  fi  ♦  B*c2 


where  gQ  is  void  geometry  in  the  case  of  non-interacting  voids.  Then  the  DS  for  a  porous  material 
assumes  the  form 


^--C*A>-v2)  (a) 

^--S*B>-v2)  (b)  (19) 

Ao-Ao(C’-80]  Bo-B5p*-^0)  <c> 

Now  multiply  both  sides  of  (19)  by  v2  and  perform  the  limit  of  flattening  out  the  voids  into  cracks. 
Then  for  example 


lim 


,  dc: 

2  dv„ 


:*lim(A*v2  ) 


Since  A*  is  a  special  case  of  Aq  it  certainly  satisfies  (14)  and  therefore  (19c)  also  satisfies  (14). 
Recalling  (17c)  it  follows  that 


lim  (A*  v2 )  -  R* 


with  a  similar  reasoning  for  B  we  can  finally  write  the  modification  of  (19)  to  the  case  of  a 


cracked  solid  in  the  form 


lim  v,  -  -  cV 

2  dv2 


tov2  %  -SV-R* 


R’-R,(S*«c) 


v2 


C*  -  C(1) 

v2-0 


s*  -  s(,) 

v2  «  0 


The  limit  on  the  left  side  cannot  be  carried  out  in  a  general  sense  and  depends  on  the  crack 
geometry  as  will  be  seen  further  below. 

We  specialize  (20)  to  the  important  case  when  the  cracked  material  is  statistically  isotropic.  In 
this  event  all  tensors  appearing  in  (20)  must  be  isotropic  and  we  represent  them  in  the  convenient 
notation.  Hill  (1965),  Walpole  (1981) 


C*  -  3K*!1  +  2G*!2 


5’  "  (3K*  '  2G*  ] 

?'-(3V2P2] 


(3K*  ,2G* ) 


& 

s 


ft 

5? 


S'  vSSSSS  sss 


m At  »«»  m  ;»»«■  *«»  m  »»»:»..< 


K!WTOJWV>JWAn.  wwirv.  irv,  i ru  irv.  »r*f  ar\*  wx*  wv»w. «  v»  >w*  *■»-»•' v^»-v»  c 


where 
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fl  j  r  c 

ijW  yijkt 


I1 11. 11 


'fjtt -  i  [Vji'Vjk  - !  Vu] 


I2!2  -I2  I1 12  -  I2!1 


(22) 


The  great  advantage  of  this  notation  is  that  a  product  such  as  in  the  right  side  of  (20)  assumes  the 
simple  form 


CV  -  pK*P*.  4G*P*  j 


Introducing  these  results  into  (20a)  we  have 


lim 


lim 


, 

2  dv„ 


,  d X 

2  dv„ 


-  -  3KV 


-  .g*p; 


(23) 


with  analogous  results  for  (20b).  Obviously,  in  the  case  of  small  crack  density  the  material  will  also 
be  statistically  isotropic.  Therefore,  the  tensors  (16c)  admit  the  representation  (21)  and  thus 


vlK^c) 


(24) 


3 
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with  a  similar  result  for  Q  . 


APPLICATIONS 


We  consider  first  the  case  of  randomly  oriented  elliptical  cracks  which  have  different  sizes  but 
identical  ratio  of  minor  to  major  axes  a/b.  The  small  crack  density  result  is  implicitly  given  in 
Budiansky  and  O'Connell  (1976)  and  is  here  put  into  the  form 


K,  -  K(1  -  /to) 


G(1  -  /icr) 


(25) 


where  K  and  G  are  matrix  properties. 


a  » 


xEab2 

V 


(26) 


is  the  crack  density  parameter  (CDP)  which  enters  from  evaluation  of  (12)  and  the  sum  indicates 
summation  over  all  cracks.  Furthermore 


*  -  x(b/a,«/) 


8(»tV2) 

9(l-2v)E(k) 


ft  -  b/z,v) 


!$  Jb±.  1  +  3  f  1 

45  E(k)  4  [  Upu  Upu  J 


(27) 


where 
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jkjK(k)-  E(k)|/k2E(k) 

P  ■  /kj[E(k)  -  K(k)]\ /k2  E(k) 


k2  -  1  -  b2/a2 


k2-b2/a2 


(28) 


Here  v  is  the  matrix  Poisson's  ratio  and  K(k)  and  E(k)  are  the  complete  elliptic  integrals  of  first  and 
second  kind,  respectively.  Note  that  the  bias  towards  b  in  (26)  is  due  to  the  choice  of  representa¬ 
tion  in  terms  of  the  classical  elliptic  integrals. 

The  equivalence  of  (16a,c)  and  (25),  utilizing  the  isotropy  of  all  tensors,  leads  to  the  identifica¬ 
tion 


1  "  3  ** 


Therefore  from  (24) 


* 


\*a 


1 

2 


M  <* 


k*  m  k  (b/aV )  ft*  -  n  (b/a ,u* ) 


(29) 


(30) 


where 


•  m  3K*  -  2G* 
"  2(3K*  +  G*) 


(31) 


is  the  effective  Poisson's  ratio.  It  remains  to  evaluate  the  limits  on  the  left  side  of  (23).  Assume 
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that  the  elliptical  cracks  are  obtained  by  flattening  ellipsoidal  voids  which  have  the  three  axes  a,  b, 
pb  where  p  can  be  made  indefinitely  small.  Then  the  volume  fraction  of  voids  is 


I 


Vj  -  ®sa_.p« 


(32) 


where  the  last  equality  is  due  to  (26).  Remembering  that  dv^  is  produced  by  adding  voids  it  fol¬ 
lows  that 


dv^  -  pda 


(33) 


Inserting  (32,33)  into  (23),  p  cancels.  Introducing  (30)  we  obtain 


dK* 

da 


dG* 

da 


-kV 


n*  * 

G  (i 


(34) 


which  in  view  of  (31)  are  two  simultaneous  differential  equations  for  K*  and  G*.  with  initial  con¬ 
ditions 


a-0 


which  can  be  easily  integrated  numerically. 


a*0 


In  the  special  case  of  penny  shaped  cracks,  a-b,  these  equations  reduce  to 


(35) 


ft 


V-i 


I 


.  Ji.  II.  I'*  »'•  *' ii  a,  h  , 


I£k.#  1^- 

6  l.9u 


32g*  (1-/XS-/) 
45  U  ?  • 


dE*  16  p*  (!->/  X10-3i/* 
da  45  •j.w* 


a  -  *£a  /V 


The  first  three  equations  are  interdependent  and  any  of  them  follows  from  the  other  two.  These 
equations  have  already  been  stated  by  Salganik  (1973)  on  the  basis  of  the  Roscoe  (1952)  interpreta¬ 
tion  of  the  DS  according  to  which  cracks  added  are  always  by  an  order  magnitude  larger  than  pre¬ 
vious  ones.  The  equations  can  be  integrated  in  closed  form  as  is  shown  in  appendix  II,  and  the  res¬ 
ults  are 


„  _  5  ,n  v  15  .  l-i/*  45  .  \*v* 

a«;en-5-+2T®  -  ♦  TVS'  U1  - 

8  u*  64  l-i/  128  l*i/ 


*  S'"  SXT 


l*i/  E 


These  expressions  have  also  been  given  by  Zimmerman  (1985). 

Fig.  1  shows  the  variation  of  G*/G  and  E*/E  with  crack  density  on  the  basis  of  (36).  Also 
shown  are  the  SCS  predictions  for  this  case.  The  latter  illustrate  the  difficulties  mentioned  in  the 
introduction  and  it  is  seen  that  in  contrast  the  DS  results  vanish  asymptotically  with  increasing  crack 
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density. 

The  calculations  just  performed  reveal  some  noteworthy  aspects.  The  first  one  is  that  in  deriv¬ 
ing  (34)  the  CDP  a  cancelled  on  both  sides  of  the  equation.  It  is  believed  that  this  is  peculiar  to 
gHiprieal  cracks  and  stems  from  the  strain  uniformity  in  the  ellipsoidal  inclusion  problem.  For 
cracks  of  more  general  shape  it  cannot  be  expected  that  the  limit  on  the  left  side  and  the  mathemat¬ 
ics  involved  in  the  right  side  of  (23)  should  result  in  the  same  geometrical  parameter. 

The  second  observation  is  that  the  DS  as  interpreted  here  and  by  Roscoe  (1932)  give  the  same 
results  for  elliptical  cracks  with  constant  ratio  a/b  but  not  for  cracks  with  different  a/b.  Therefore 
the  two  approaches  are  not  necessarily  equivalent. 

The  second  example  of  DS  application  is  concerned  with  the  stiffness  reduction  of  a  plane  orth¬ 
otropic  thin  layer  due  to  a  distribution  of  parallel  line  cracks,  fig.  2.  This  may  be  regarded  as  a 
possible  model  for  a  cracked  ply  inside  a  laminate  consisting  of  unidirectional  fiber  composite  plies 
and  it  is  therefore  of  interest  to  consider  the  case  of  constant  plane  stress  j »  •  a\2 '  Consider 

first  the  case  of  a  single  crack  of  length  2a  within  an  orthotropic  ply,  fig.  3.  Stress-strain  relations 
in  plane  stress  are 


(38) 


€  .  11 

12  2G12 


It  is  easily  realized  that  application  of  <7j  j  by  itself  does  not  produce  any  COD.  When  only  is 
applied  the  COD  (Uj  ]  vanishes  because  of  symmetry  and  when  only  is  applied  the  COD  [u^  ] 
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vanishes  because  of  antisymmetry.  The  surviving  COD  are  well  known;  see  e.g.  Sih  and  Liebowitz 
(1970),  and  may  be  written  in  the  form 


I“l  1  “  *^12  n  J*2  •  *2 

[“2 1 "  °22  A  J®2-  *2 

(39) 

2>/5  r  1/2  1I/2 

*  -  ♦V2G.2-*'>2] 

A  -  (E,)/(E2),/2  7 

Next  we  evaluate  as  given  by  (12)  for  the  case  of  small  crack  density,  in  which  case  all 
cracks  have  the  COD  (39).  Noting  that  for  the  present  case  ij-1,2  and  n^  -0,  n2-l,  it  follows 
easily  that 


'11 


„  .  Aq 

722  2  22 


712  "  4**12 


(40) 


where  a  is  the  CDP  given  by 


(41) 
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a  is  half  crack  length  and  A  is  the  area  of  the  ply  considered  of  unit  thickness.  Identifying  (40) 
with  (13b)  we  realize  that  the  only  surviving  components  of  the  tensor  Rs  of  (16)  are 


nS  AQ 

2222  2 


Rs  - 
1212  S 


(42) 


The  effective  stress-strain  relations  of  the  cracked  ply  are 


*  11  "Sll!l  *11  +S1122ff2 2 


*22  “SU22aH  +S2222  *  22 


(43) 


*  12  “2S12I2*12 


where 


i 

i 


(44) 


To  set  up  the  DS  we  employ  (20b,c).  We  consider  the  line  cracks  as  limiting  case  of  elliptical 
cracks  with  major  axis  a  and  minor  axis  pa.  Then 


v 


2 


pxEa^ 

V  " 


pa 


(45) 


dv2  -  pd 


From  (20c)  and  (42)  the  only  surviving  components  of  R~  are 


7777  2 


*  • 

R  - 

1212  8 


where  the  asterisk  implies  that  all  properties  in  rj  and  n  of  (39)  are  to  be  replaced  by  effective  pro¬ 
perties.  It  now  follows  from  (20b),  (43-46)  that 


2222 

da 


Eqs.  (47a)  imply  that  these  properties  are  not  function  of  crack  density  and  therefore  they  are  equal 
to  the  properties  of  the  uncracked  material.  Thus 


VE1 


* 

v  m  v. 
12  1 


This  is  of  course  an  exact  result  which  is  easily  deduced  from  first  principles.  It  follows  that 


•  2^2  T  *1/2  •  V 

7  -  e7[(VV  +e./2G.2-^ 


'A  A’.V  s'j.w ' n' ^  -.'IV.'v' \ ■‘•'S'InV *  ' 


LTWWWl'V 


’.VWWATV  V  w  •- >  i sr  V-  v 


and  therefore  from  (44,47) 
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A  -  (Ej/Ej)  9 


dE 

—  S2 

da 


.3/2 


1/2 


1/2 


2G 


12 


12 


1/2 


dG 


12 


da 


-- V2 


,•  2 
*12 


1/2 


2G 


12 


1 1/2 


12 


(49) 


(50) 


Eqs.  (50)  are  easily  integrated  in  closed  form.  If  the  equations  are  divided  one  by  another  we 

•  * 

obtain  a  simple  differential  relation  between  E_  and  G  _  whose  integral  is 

2  12 


i  i 

o;2 ' 


(51) 


Expressing  from  (51)  and  substituting  into  the  second  (50)  we  obtain  a  simple  differential  equa¬ 
tion  with  integral 


r*.i 

1/2 

.  1 

fE'l 

1/2 *  E-  1 

K 

E2 

2°  12  12 

E2 

2°  ,2  ^ 

1/2 


i/y/2 


(52) 


It  is  easily  seen  that  (S3)  imply  the  relation 


r  • 

1/2 

1  _ 

E2 

1  1 

°12 

lE>] 

<  "  E2 

* 

(54) 


•  * 

and  substituting  for  from  (54)  into  (53)  we  obtain  an  equation  for  E^. 

Fig.  4  shows  plots  of  e'/E^  and  as  functions  of  CDP  for  a  graphite/epoxy  ply, 

according  to  the  DS  and  the  SCS.  Note  that  in  this  case,  unlike  the  case  of  penny  shaped  cracks, 
there  is  no  dramatic  difference  between  the  two  predictions.  Note  also  that  shear  modulus  reduction 
is  much  smaller  than  transverse  Young's  modulus  reduction.  This  is  apparently  due  to  the  large 
stiffness  in  fiber  direction  which  inhibits  the  shear  COD  which  is  also  in  fiber  direction 


CT\PTilr  HJWI7  ,TB  nm  rwui 
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CONCLUSION 

A  differential  scheme  (DS)  approximation  to  obtain  the  elastic  properties  of  solids  containing 
distributions  of  many  cracks  has  been  established  in  a  general  sense.  The  necessary  information  for 
application  of  the  method  is  the  COD  of  a  typical  single  member  crack  of  the  distribution  when 
isolated  in  infinite  homegeneous  medium  with  elastic  symmetry  identical  to  the  macro  symmetry  of 
the  cracked  material. 

Specific  results  have  been  derived  for  randomly  oriented  elliptical  and  penny  shaped  cracks  in 
elastic  matrix  and  for  aligned  plane  cracks  in  orthotropic  sheet  There  is  a  drastic  difference 
between  the  DS  and  SCS  results  for  penny  shaped  cracks  and  much  milder  difference  for  aligned 
plane  cracks. 

The  method  can  be  readily  applied  to  many  other  cases  of  interest 
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APPENDIX  I 


SELF-CONSISTENCY  OF  DIFFERENTIAL  SCHEME 


It  is  required  to  prove  that  C*  and  S*  as  predicted  by  (5,6)  are  mutually  reciprocal.  Thus 


cV  - 1 


(i-i) 


and  therefore. 


S*  +  C*  m  o 
dv2  dv2 


Substituting  from  (5,6)  into  the  second  of  (1-1)  we  have 


<C(2)  -  c')  A*S-  ♦  C*(S(2)  -  sW 


(Sv  /  -  S  )B  -  0 


a-2) 


Recall  the  definition  (4)  of  A*  and  analogous  definition  of  B*  as  small  concentration  tensors  in 
which  matrix  properties  have  been  replaced  by  effective  properties.  The  left  eqs.  (3)  must  be  mut¬ 
ually  reciprocal.  Carrying  out  the  multiplication  to  order  c2  we  have 


(C<2)  -  C(1))  A  S(1)  +  C(I)(S(2) 
$ 


)Bs-0 


(1-3) 


Replacing  C^  and  S^  by  C*  and  S*  in  (1-3)  proves  (1-2)  and  thus  the  second  of  (I- 1 ).  This 
implies  that  the  product  C*S*  is  a  constant  tensor  which  must  however  be  equal  to  I  because  of  the 
known  reciprocity  of  C^  and  S^  in  the  initial  conditions  of  (5,6). 

This  completes  the  self-consistency  proof. 


APPENDIX  II 


SOLUTION  OF  THE  DS  EQUATIONS  FOR  RANDOMLY  ORIENTED  PENNY  SHAPED 
CRACKS 


Write  the  second  and  third  eqs.  (36)  in  the  form 


*  F(v*  )da 

Ct 

^--H(»/*)da 

G 


(0-1) 


By  definition 


E*  -  2(1  +  j/*)G* 


(H-2) 


Combination  of  the  differential  of  (II-2)  with  (II- 1)  yields 


da 


_ d/ _ 

(l+i/*XF(i/*)  -  H(i/* )] 


(H-3) 


and  thus 


f  di * _ 

J*  (l+i/* H(i/*)] 


(H-4) 


Introducing  (II- 3)  into  (II- 1)  with  the  definitions 


v 


l 

73 

.'■I 
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e  -  E*/E 


g  -  G  /G 


we  have 


toe  -  r  _ ElL^iL _ 

)v  (1+/XF(^)  -  H(i/*)J 

to  g  -  T  _ HfeV* _ 

J„  (l+y*XF(y*)-H(y*)) 


Introducing  F(j/*)  and  H(i/*)  from  (36)  the  results  (37)  follow  by  elementary  integration.  There  is 
of  course  no  difficulty  to  employ  (11-4,6)  to  analyze  the  case  of  elliptical  cracks  in  accordance  with 
(25)  and  (27,28).  In  this  case  it  is  best  to  carry  out  the  integrations  numerically. 
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APPENDIX  II 


SOLUTION  OF  THE  DS  EQUATIONS  FOR  RANDOMLY  ORIENTED  PENNY  SHAPED 
CRACKS 

Write  the  second  and  third  eqs.  (36)  in  the  form 


~~  -  T(u*  )da 


H(/)do 

G 


By  definition 


E*  -  2(1  +  v*Yj* 


Combination  of  the  differential  of  (13-2)  with  (II- 1)  yields 


da 


d v 


(1+/XF(^)  -  H(y* )] 


and  thus 


r 

iu  (1+J/* XF(v* )  -  H(1/* 


)] 


Introducing  (II-3)  into  (II- 1)  with  the  definitions 


(H-l) 


(n-2) 


(n-3) 


(n-4) 


A 


S.1 


a 


v" 
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